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ABSTRACT 



We present results from six years of recalibrated and new spectropolarimet- 
ric data taken with the University of Wisconsin's Half- Wave Spectropolarimeter 
(HPOL) and six years of new data taken with the photoelastic modulating po- 
larimeter (PEMP) at the Flower and Cook Observatory. Combining these data 
with polarimetric data from the literature allows us to characterize the intrinsic 
BVRI polarized light curves. A repeatable discrepancy of 0.245 days (approxi- 
mately 6 hours) between the secondary minima in the total light curve and the 
polarization curve in the V band, with similar behavior in the other bands, may 
represent the first direct evidence for an accretion hot spot on the disk edge. 

Subject headings: (Stars:) binaries: eclipsing - stars: individual ((3 Lyrae) - 
techniques: polarimetric - accretion, accretion disks 



1. INTRODUCTION 



Beta Lyrae A (also known as HD 174638, HR 7106 and ADS 11745A; hereafter "/3 
Lyr") is a bright, well-studied semi-detached eclipsing binary star system. The primary star 
is a B6-B8 II, giant star ("loser") with a mass of 3 M that is transferring matter to its 
main-sequence BO. 5 V, 12.5 Mm, companion ("gainer") at about 10 -5 M n yr" 1 via Roche 



lobe overflow ( 
a thick accreti< 


Hubenv & Plavec 


1991; 


Harmanec & Scholz 


993). This process has created 


m disk that obscures the gainer ( 


Huang 


1963; 


Wilson 


1974; 


Hubenv & Plavec 



199ll ; ISkulskiilll992l ). A bipolar flow or jet has a lso been detected in the system through in- 



terferometric and spectropolarimetric methods (IHarmanec et al.lll996l ; iHoffman et al.lll998 
hereafter HNF). The s ystem's mass ratio, q, has been placed between 4.2 and 6 with an incli- 
nation angle, i, of 85° fl Wilson! 1 1974T). Other studies have sugg ested i = 83° with q = 5.6 and 



i = 80° with q = 4.28 (IHubenv fc Plavedll99ll : ISkulskii 
of the orbital inclination place its value at % 



1992). More recent determinations 



86° flLinnell et al.lll998l : iLinnelll 120001 ). These 
large mass ratios are evidence for mass reversal in the system's history. The disk's ability to 
obscure the gainer is due to the nearly edge-on inclination angle of the system. 

The system has a well-establish ed orbital period of 12.9 days that increases at a rate of 
19 s yr -1 (IHarmanec fc Scholzlll993l ). Recent interferometric observations have produced the 
first images of the system, which show the loser and the disk as separate objects and confirm 
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the orien tation of the system axis, near 254°, previously inferred from HNF's polarimetric 



analysis ( IHarmanec et al.lll996l ; IZhao et al.ll2008l ; ISchmitt et al.ll2009l ). Understanding how 



mass moves between and around the stars and leaves the system is imperative to understand- 
ing the evolutionary future of j3 Lyr. However, interferometric techniques have yet to resolve 
the mass stream or bipolar outflows. We have used spectropolarimetry to study the system 
and begin to unlock the evolutionary clues contained within the circumstellar material. 

Light scattering from electrons in the highly ionized circumstellar material in (3 Lyr 
produces a variable phase-dependent polarization. Since electron scattering preserves infor- 
mation about the orientation of the scattering region, analyzing polarimetric behavior as a 
function of wavelength allows us to determine from where in the system different spectral 
features have arisen. In this way, spectropolarimetric observations of /3 Lyr can be used to 
infer the geometrical properties of the scattering material in the system. 

Optical polarimetry was used to study Lyr as earl y as 1934, but it was not know n 
until 1963 that the system exhibited variable polarization (jOhmanlll934l ; IShakhovskoilll963l ). 
Appenzeller and Hiltner (1967; hereafter AH) were the first to publish interstellar polar- 
ization (ISP) corrected broadband UBV polarization curves of /3 Lyr. More recently, HNF 
published ISP-corrected polarized light curves in the V band and Ha and He I A5876 emis- 
sion lines using a subset of the data we present here. HNF used the position angles of the 
polarized UV continuum and the hydrogen Balmer emission lines to confirm that a bipolar 
outflow exists in the /3 Lyr system after their discovery by Harmanec et al. (1996). HNF also 
interpreted the average position angle of the visible polarized light (164°) to be the physi- 
cal axis of the binary system, an interpretation which was borne out by the interferometric 
images presented by Zhao et al (2008) and Schmitt et al. (2009). 

In this paper we present new BVRI and He I A5876 polarization curves and polarized 
light curves of /3 Lyr. The details of our spectropolarimetric observations and our interstellar 
polarization corrections are in Section 2. Section 3 presents and displays our observational 
results. We analyze our findings in Section 4 and summarize conclusions in Section 5. 



2. OBSERVATIONS 

This study compiles data from three distinct data sets. The first consists of 69 optical 
spectropolarimetric observations of /3 Lyr taken over 6 years with the University of Wiscon- 
sin's Half- Wave Spectropolarimeter (HPOL) at the 0.9 m telescope at Pine Bluff Observatory 
(PBO); the second data set comprises 6 years of broadband optical polarimetric data ob- 
tained at the Flower and Cook Observatory; and the third is 3 years of archived broadband 
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optical polarimetric data from AH taken with the 24 inch rotatable telescope at the Yerkes 
Observatory. To calculate the phase for each observation, we used the ephemeris 

T pri = HJD 2, 408, 247.966 + 12.91378E + 3.87196 x l(r 6 £ 2 

where E is th e total number of orbits since the primary eclipse that occurred at HJD 
2,408,247.966 (IHarmanec fc Scholzl Il993l ). This is the same ephemeris used by HNF; it 
does not significantly differ from the more recent ephemeris presented by Ak et al. (2007). 



2.1. HPOL Data Set 



The first 14 HPOL observations, obtained between 1992 September and 1994 November, 
used a dual Reticon array detector with a wavelength range of 3200-7600 A and a resolution of 
15 A (see Wolff, Nordsieck, & Nook 1996 for further instrument information). The remaining 
55 observations, taken between 1995 March and 1998 September, used a CCD-based system. 
This extended the wavelength r ange, 3200-10500 A, and increased the resolution to 7.5 A 
below 6000 A and 10 A above (INordsieck fc Harris! Il996l ). The first 29 observations were 
previously published in HNF; they have undergone recalibration for use in this study. 

Table [1] lists the orbital phases along with civil and heliocentric Julian dates for the 
midpoint of each HPOL observation. Each /3 Lyr observation covers the full spectral range, 
with the exception of the four nights indicated in Table [TJ Two observations, 1995 May 27 
and 1995 August 14, used only the red grating (6,000-10500 A) of the CCD system while the 
other two, 1997 May 17 and 1997 May 26, used only the blue grating (3200-6000 A). Each 
individual observation typically lasted between 45 minutes and an hour (approximately 0.03 
to 0.04 days) when both gratings were used. 

We used 11 HPOL observations of /3 Lyr B taken between 1995 May and 1999 November 
to obtain an ISP esti mate. Beta Lyr B is a member of the same association as ft Lyr A and 
is located 45" away ( lAbt et al.lll962l ). All of these observations were made using HPOL's 
CCD-based system. Table [T] also lists the civil and heliocentric Julian dates that correspond 
to the midpoints of each observation of f3 Lyr B and indicates which grating(s) were used 
during the observations. The first two observations, 1995 May 21 and 1996 July 3, were 
previously published in HNF and have undergone recalibration for this study. We reduced 
all of the HPOL observations using the REDUCE software package (described by Wolff et 
al. 1996). 

We estimated the ISP by fitting a modified Serkowski law curve to the error-weighted 
mean of the 11 observations of /3 Lyr B (jSerkowskil Il965l ; IWilking et al.lll982l ). The param- 
eters for our ISP estimate are P max = 0.422% ± 0.005%, A max = 4149 A ± 80 A, K = 0.699 
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and PA = 151.16° ± 0.36°. We subtracted this ISP estimate from the HPOL Lyr data. 
This new ISP correction has significantly improved uncertainties over previous estimates; 
it is consistent with the HNF estimate (P max = 0.419% ± 0.013%, A max = 4605 A ± 260 
A, and P.A. = 151.0° ± 0.9°), which was determined using only the 1995 May 21 HPOL 
observation. AH found a similar estimate, P max = 0.42% ± 0.04% and PA = 153.2° ± 3° 
by taking the weighted mean of the observed polarization of the associated stars (3 Lyr B, 
E, and F. 



2.2. FCO Data Set 



Our second data set is made up of 19 B band, 88 V band and 17 R band observations 
obtained at the Flower and Cook Observato ry between 1987 and 1992 with the PEMP in- 
strument (IHolensteinl Il99ll ; lElias et al.lll996l ). The length of each observation was between 
20 and 25 minutes (approximately 0.01 days). The phases along with the civil and heliocen- 
tric Julian dates for each observation are listed in Table [21 We have no observations of (3 
Lyr B taken with the same instrument as this data set. Therefore, we used the Serkowski 
fit to the HPOL (3 Lyr B observations (see Section 2.1) to calculate the ISP contributions 
at the central wavelengths of the BVR bands and subtracted these estimated values from 
the observations in this data set. We list these data and ISP subtracted data in Tables [3] 
through [51 



2.3. AH Data Set 

We also used archival BV polarization data taken between 1964 and 1966, originally 
published in AH. This data set consists of 37 B band and 127 V band observations, the details 
of which can be found in AH (and references therein). The DC two- chan nel polarimeter wa s 



rotated 30° between 12 separate 20 second exposures of (3 Lyr and the sky (lAppenzellerlll965l ) . 
Therefore, the total integration time for both (3 Lyr and the sky was 4 minutes (approximately 
0.003 days). We converted these data from polarization magnitudes to percent polarization 
and converted their Julian dates to heliocentric Julian dates for use in this study. For 
consistency between the HPOL and AH data sets, we did not use AH's published ISP 
corrected data because it included two stars ((3 Lyr E and F) that are not included in the 
HPOL ISP estimate. Instead, we subtracted only the AH BV (3 Lyr B observations from 
their non-ISP corrected (3 Lyr data in each respective band. 



-6 - 



3. RESULTS 

3.1. Broadband Polarimetry 

To examine the behavior of the continuum polarization with orbital phase, we ap- 
plied synthetic BVRI Johnson- Cousins band filters (described by Bessell 1990) to the ISP- 
corrected HPOL data. The filter routine produces broadband values and associated internal 
errors for each observation; however, we must still take into account systematic variations 
in the instrumental polarization between nights. Systematic errors for HPOL at PBO were 
evaluated by periodically analyzing observations of unpolarized standard stars. Tables [6] 
through [9] list broadband polarization values and internal errors determined by the filter 
routine along with the systematic errors. In the case of the Reticon data, the systematic 
errors are less well determined; based on our previous experience with these data, we have 
estimated the uncertainties in the Stokes parameters for the Reticon data to be 0.02% in all 
bands. Figures [1] through H] display these data graphically, using the larger of the internal 
and systematic errors for each observation, along with the FCO and AH data. 

The position angles for the BVRI bands remain relatively constant with orbital phase 
except near secondary eclipse, where the position angle values appear to rotate away from 
the mean value. The bottom panels in Figures [1] through H] show this behavior in each band. 
In the standard picture of the system, the polarization is produced by scattering from the 
accretion disk edge; the average position angle at these wavelengths should therefore provide 
us with an estimate of the orientation of the axis of the disk and thus of the system as 
a whole. To calculate the average position angle in each band, we excluded the secondary 
eclipse points (between phases 0.425 and 0.575) because they do not follow the near-constant 
trend displayed at other phases and performed a linear, error-weighted, least-squares fit to 
the remaining data in Q-U space. We then used the slopes of the fitted lines to determine 
the position angles listed in Table [TU1 We also calculated a position angle for Balmer jump 
index (the vector difference between the polarization above and below the Balmer jump). 
As discussed in HNF, the broadband polarization in (3 Lyr undergoes a 90° position angle 
rotation across the Balmer Jump. Thus, this vector difference defines the system axis in 
Q-U space. Because the Balmer jump index is independent of the ISP, it provides us with 
an independent estimate of the orientation of the system axis. 

We find the weighted mean position angle of the HPOL Reticon and CCD bands and 
the Balmer jump index to be 164.6° ± 0.22°. For each band, the CCD and Reticon position 
angles do not agree within uncertainties, with the Reticon data yielding larger position 
angles in all bands. This result is due to poor sampling in Q-U space; the Reticon data only 
consist of 14 data points in each band while there are 55 observations in the CCD V and 
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R bands and 53 observations in the CCD B and / bands. This skews the linear fit since 
the full range of possible observable Q-U values is not well covered by the Reticon data. 
However, the larger systematic uncertainties we adopt for the Reticon data result in these 
points carrying a lower weight in the fit; thus, we are confident that our weighted mean is 
a fair representation of the true system axis. HNF found the mean V band position angle 
to be 163.8° ± 0.15°; while our estimate is not formally consistent with HNF's, it shows 
the broad band polarization behavior of the system is the same in all optical bands. Our 
polarization position angle implies a system position angle, defined by the position angle of 
the disk axis, of 253.8° ±0.15° on the sky. As expected for polarization by electron scattering, 
the position angle of the polarized light is perpendicular to the position angles describing 
the system orientation given by Zhao et al. (2008) and Schmitt et al. (2009). Zhao et al. 
(2008) estimated the position angle of the system's ascending node as 253.22° ± 1.97° and 
251.87° ± 1.83° using two different image reconstruction techniques on their interferometric 
data and 254.39° ± 0.83° using a model of the system, while Schmitt et al. (2009) estimate 
249.0° ±4.0°. 

We rotated all of the HPOL, AH and FCO data to the average position angle of 164°. 
This orients our data with respect to the intrinsic polarization axis of the system. After this 
rotation, % U averages to zero in each band and the polarization varies significantly only in 
the %Q direction. In the rest of this paper we present the projected Stokes parameter %Q p 
resulting from this rotation. The use of this quantity is beneficial because it can be positive 
or negative, whereas %P is always positive. Data points that have a position angle near 
164° will have a positive % Q p value while points with a position angle perpendicular to this 
(near 74°) will have a negative %Q p . Hereafter we display only %Q p because %U P values 
scatter around zero. Because the rotation is a simple trigonometric calculation, we present 
in Tables [3] through [9] the unrotated % Q and ISP subtracted FCO and HPOL data only. 

The middle panels in Figures [T] through H] show the %Q p curves for the BVRI bands after 
rotation. We used the program PERIOD04 to perform a Fourier fit to the data for each band 



flLenz fc Bregeri 120051 ). The PERIOD04 fitting formula is y = Z + M sin^vr^t + <&)) 
where n is the number of sine terms in the fit, Z is the zero point, A is the amplitude, Q 
is the frequency, and (f> is the phase. The results of the fits are displayed as solid curves in 
each figure and their parameters are given in Table [TTJ The V band, and to a lesser extent 
the B band, Fourier fits deviate from the data at phase 0.9 (see Figure [2]). This discrepancy 
disappears in the V band if we include three frequency terms in the Fourier fit. However, we 
are not as confident in the third frequency as we are in the first two because the PERIOD04 
fitting program produces a reasonable third term for the V band only. Therefore, in Table 
[TT1 we report parameters for only the first two terms, but we display both the two-term and 
three-term fits in Figure [2j These fits provide the first quantitative representations of the 
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polarization variations in the (3 Lyr system. 

The data in the BVRI bands are almost always positive, indicating their position angles 
stay near 164° throughout the orbital period. Each of the bands displays an increase in %Q p 
at primary eclipse and two other increases near the quadrature phases (0.25 and 0.75). The 
height difference between the polarization bumps at the quadrature phases noted by HNF 
disappears now that more data are added, but we note that the B and V phase 0.25 bump 
have a higher dispersion around the average %Q p value than does the 0.75 phase bump. 
The R and / band show the opposite behavior; future observations will be able to tell us 
whether this is due to their poor phase coverage or whether it indicates that the R and / 
bands are probing a different region of the disk than the B and V bands. We calculated the 
variance of the two quadrature bumps between phases 0.25 and 0.35, and 0.65 and 0.75 to 
formally show this. The first quadrature bump has BVRI variances of 0.081 ± 0.006, 0.110 
± 0.005, 0.040 ± 0.002, and 0.016 ± 0.002 respectively, while the second quadrature bump 
has variances of 0.048 ± 0.005, 0.027 ± 0.002, 0.281 ± 0.005, and 0.069 ± 0.004. 

The R and / bands also produce a lower polarization signal than the B and V bands. 
This change in polarization behavior with wavelength could indicate that scattering mecha- 
nisms other than electron scattering are present in the system. However, more observations 
in the R and / bands are needed to rule out the possibility that the low signal is due to a 
lack of phase coverage in these bands. 

Figures [T] through H] also show that in each band, the polarization curve has a fitted 
minimum in polarization that occurs just prior to the secondary eclipse in total light. This 
minimum is accompanied by a rotation in the position angle of the polarized light away from 
its average value. Table [121 gives the phases of the %Q P Fourier fit minimum in all bands. 
The offset between secondary eclipse in polarized light and total light is a new result, seen 
here for the first time due to the improved phase coverage in these data. We discuss the 
implications of this phenomenon in Section 4. 

At primary eclipse, there are hints of similar behavior: the polarization maxima in the B 
and V bands occur slightly before phase 0.0, and all bands show a deviation from the mean 
position angle at and just after phase 0.0. However, we consider the primary eclipse features 
to be less significant than the ones at secondary eclipse for the following reasons. The three 
data points showing noticeable deviations from the mean position angle all occurred on the 
same night, 1997 August 25, which suggests that this effect may be due to a non-periodic 
process intrinsic to the Lyr system, or to a change in observing conditions that affected 
that night's data. If the same structure is responsible for the phenomena at both eclipses, 
we expect that a similar position angle scatter should exist in observations from the same 
orbit of the system whose phases are between primary and secondary eclipse, when such a 
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structure should be most visible to the observer. However, neither the data from 1997 August 
26, nor 1997 August 30 show such behavior. Additionally, the polarization maximum in the 
R band occurs at phase 0.0 within the uncertainties, and in the I band the maximum occurs 
just after phase 0.0. If we include the third term in the V band Fourier fit, the primary 
polarization maximum occurs at phase 0.0 within uncertainties. Because the eclipse behavior 
is not consistent between bands and the associated position angle scatter appears to have 
occurred during only one orbit of the system, we consider it unlikely that these are due to a 
stable physical structure within the system (see Section 4 for further discussion of the effects 
seen at primary eclipse). 

Figure [5] displays the projected polarized BVRI flux light curves for the (3 Lyr system. 
To create these polarized light curves, we multiplied the fitted pol arization curves shown 
in Figures [T] through 0] by their respective Fourier fit light curves (IHarmanec et al.l Il996l ) 
normalized to maximum light. In all bands, the polarized flux remains nearly constant 
across primary eclipse due to a decrease in total light and an increase in %Q P . The secondary 
eclipse offset seen in Figures[T]through|4]persists in Figure[5j while quadrature phases display 
local maxima. The B and V band appear to be the most similar; they overlap for most 
phases, while outside of secondary eclipse and the first quadrature phase the R and I bands 
produce the lowest net projected polarized flux. We do not consider the apparent height 
differences between bands at the quadrature phases to be significant due to the scatter in 
the observational points and the lack of coverage in the R and / bands. 



3.2. Line Polarimetry 

We also took advantage of the spectropolarimetric nature of the data by studying the 
polarization behavior of the strongest optical emission lines in /3 Lyr's spectrum: Ha, H/3, 
He I A5876, He I A6678, and He I A7065. HNF hypothesized that the Ha, H/3, He I A5876 
and He I A7065 lines, which show a negative projected polarization, scatter in the bipolar 
outflow, while the He I A6678 scatters on the edge of the disk. However, HNF did not have 
enough data to construct a full polarization phase curve for the lines. Our expanded data 
set allows us to do this. However, we present only our He I A5876 results in graphical form 
because our uncertainties are relatively large due to signal-to-noise limitations. Rather than 
present the Ha, B./3, He I A6678, and He I A7065 data, we describe their general behavior 
below. Future observations will allow us to use these data to draw quantitative conclusions 
about the scattering regions that give rise to the polarization in these lines. 

In order to calculate the polarization for each emission line, we used the flux equiva- 
lent width method described by HNF, using the same line and continuum regions as far as 
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possible. We corrected the Ha and H/3 lines for underlying unpolarized absorption compo- 
nents (arising from the loser) in the same manner as HNF, using their preferred absorption 
equivalent widths of 8 ± 2 A for Ha and 6 ± 1 A for H/3. This has the following effect on 
the data. The continuum is positively polarized while the lines are negatively polarized. If 
we do not correct for unpolarized absorption, we remove too much continuum, and thus our 
resulting line polarization is too negative. With the absorption correction, the continuum 
contribution is smaller and the magnitude of the polarization is also smaller, resulting in a 
less negative %Q p . 

We do not present the HPOL Reticon line polarization values due to their large uncer- 
tainties. Figure [6] shows %Q p and position angle curves for the He I A5876 line. It has a 
negative %Q p ; thus its position angle is perpendicular to the intrinsic axis of the system. 
In addition, the polarization for the line approaches zero at both primary and secondary 
eclipses. For this to happen, the scattering region for this line must lie near enough to the 
orbital plane of the system to be occulted both by the loser and by the disk. HNF previously 
suggested that this line scatters in the bipolar outflows because its average position angles lie 
near 74°, corresponding to negative values of %Q P . The results from our extended data set 
support this interpretation and further suggest that the He I A5876 scattering region within 
the outflows must lie between the loser and the disk and have a vertical extent comparable 
to the height of the disk. 

The Ha, H/3, and He I A7065 lines also all display a negative %Q p and are likely scattered 
in the same region as the He I A5876 line. Their average position angles are listed in Table 

EB 

The He I A6678 data show a polarization behavior different from that of the other lines. 
The data are generally positively polarized and their average position angle, 138.2°, agrees 
more closely with the intrinsic axis of the system than do those of the other emission lines 
(see Table HUD. 



3.3. Period Analysis 



Besides the primary orbital period of /3 Lyr (12.9 days), analysis of light curves has 
revealed several longer periodicities. A 340-day period was detected by Peel (1997), while 
both Van Hamme et al. (1995) and Harmanec et al (1996) dete cted a 282-day pe riod. Wilson 
and V an Hamme (1999) searched polarimetry from AH, HNF. ISerkowskil (jl965[ ). and IShulov 
( 119671 ) for periodicities but did not detect anything significant. 



We performed a Lomb-Scargle power spectrum analysis to search for any periodic be- 
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havior in our polarimetric broadband (HPOL, AH, and FCO) and line data (Ha and He I 
A 5876) not associated with the 12.9 day orbital period of the system. While none of the 
previously detected longer periods were found, our analysis indicates the presence of periods 
of approximately 4.3 days in both the V and B bands with a False Alarm Probability (FAP; 
Home & Baliunas 1986) of 10~ 6 when using all three data sets. None of the other bands 
or the line polarization data appear to contain periods other than the orbital period. We 
also searched for periods within the HPOL, AH and FCO data sets individually to look 
for any transient periodic variations. We find that the 4.3 day period also appears in the 
V and B AH data set with a FAP of 1CT 6 , but it does not appear in the other two data 
sets. This period is exactly one third of the 12.9-day orbital period of the /3 Lyr system 
and results from the combination of two separate effects: the increase in polarization at the 
quadrature phases due to light scattering off the disk edge and the increase in polarization 



at pr imary eclipse due to the occultation of unpolarized light by the loser ( IHoffman et al. 
20031 . hereafter HWN). These two effects cause the %Q p curves to form a complete cosine 



curve between phases 0.0 and 0.3, a second cosine curve between phases 0.3 and 0.6 and a 
third cosine between 0.6 and 1.0 (see Figures 1 through 4). Therefore, this period does not 
provide new information about the (3 Lyr system. 

We also performed a much simpler analysis to search for signatures of the 282-day 



(Wilson 


1974; 


Ak et al. 


2007) 



/0< 



Q p , %U P , position angle, 



and percent polarization for each of the BVRI bands and the He I A5876 line did not 
reveal the 282-day period and are therefore not shown. However, the time coverage of our 
observations is very uneven and could have prevented us from detecting variations on this 
time scale. The observations in our data set were taken between 1964 and 1998, but we 
have no polarimetric observations taken during 1990, 1991, or between 1966 and 1987. The 
number of observations performed per year only adds further complications; several years 
have very few observations. We also note that the only the HPOL data set has observations 
of every band and the He I A5876 line. For these reasons, we also performed an analysis 
similar to Harmanec et al. (1996) on the %Q p data for each band and the He I A5876 line. 
Plotting the data on the 282.37 day period mitigates the lack of time coverage by "folding" 
observations onto one cycle. 

We plotted the data from selected orbital period phase bins versus their phase on the 
282.37 day period (see Harmanec et al. 1996 for a similar analysis of a large amount of 
V band total-light photometry). We chose the following five orbital phase bins for several 
reasons: the 0.0-1.0 bin allows us to use all of the data, the 0.6-0.15 bin is where Harmanec 
et al. (1996) detected the 282.37 day period the most strongly, the 0.25-0.35 and 0.65-0.75 
bins allow us to determine whether the polarization of the quadrature phases changes on a 
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282.37 day time scale, and the 0.425-0.575 bin allows us to determine whether the secondary 
eclipse has a 282.37 day periodic behavior. The 0.25-0.35, 0.425-0.575, and 0.65-0.75 orbital 
period bins are too tight to leave a useful number of observations in the R and I bands and 
the He I A5876 line. While the B and V bands are slightly less affected by the size of the 
0.25-0.35, 0.425-0.575, and 0.65-0.75 orbital phase bins, their behavior does not indicate, by 
eye or by using PERIOD04, the presence of a variation on a 282.37 day period. The size 
of the 0.6-0.15 and 0.0-1.0 bins present the best chance of detecting this period because the 
amount of data is not severely reduced. However, neither bin reveals the presence of the 
282.37 day period. We also subtracted the Fourier fits from the data in the BVRI bands 
(see Section 3.1) and searched the residuals for the 282.37 day period with the same bins 
used for the original % Q data. The analysis on the residuals produced similar results; we do 
not find evidence of a 282.37 day period in either the %Q p data or the Fourier fit subtracted 
residuals. However, we note that the Fourier fits deviate from the data at some phases (see 
Section 3.1). Plots from this analysis resemble scatter plots and are therefore not shown. 



4. DISCUSSION 

In interpreting their polarized flux curves, HNF proposed two different possibilities for 
the origin and scattering location of the visible light. In their "disk-disk" case, this light 
arises from within the disk and scatters from the disk edge; in the "loser-lobe" case, the V- 
band continuum light arises from the loser and scatters from material between the loser and 
the disk. In this analysis, HNF implicitly assumed that all features of the visible polarized 
flux curve are due to the same origin and scatterer. However, modeling work by HWN and 
subsequent modeling results (Hoffman et al., in prep.) have shown that the scattered light can 
originate both from the loser and from the disk in differing proportions over the binary cycle. 
These newer results suggest the following interpretations of our BVRI polarization curves. 
The net increase in %Q p at primary eclipse (Figures [QII]) is the result of the unpolarized light 
from the primary star being blocked by the disk material at phase 0.0. HNF interpreted the 
increase in %Q p at the quadrature phases as arising in one of two ways: light originating 
from within the disk and scattering from the disk edge, or light originating from the loser 
and scattering from material between the stars. We propose, based on recent modeling work 
by Hoffman, that these "quadrature bumps" form simply by loser light scattering from the 
disk edge. The minimum at secondary eclipse occurs because the unpolarized primary star 
blocks light scattered in the secondary component. 

Near secondary eclipse in all four broadband % Q curves, the minimum in polarization 
precedes the minimum in total light; the phases for the polarization minimum in the BVRI 



13 



bands are listed in Table [12j In each band, this minimum corresponds to a rotation in 
position angle away from the average value; the phase ranges for this rotation are also listed 
in Table H2J In the basic star-star-disk model for the system, there is no mechanism to 
produce this disparity. If the loser is an unpolarized source, as indicated by the absence 
of a primary eclipse in the polarized flux curves (Figure 5; HNF), then the polarization 
minimum produced by its transit across the disk should be centered at flux minimum (phase 
0.5; HWN). Thus, to explain this offset, we need to invoke another system component. Since 
f3 Lyr is a mass transfer system, it most likely contains a mass stream connecting the loser 
and the disk as well as a 'hot sp ot' where the mass stream from the loser interacts with 



the disk edge ( ILubow fe Shulll975l . see also the geometries proposed by HNF). Some studies 
(for example, Bisikalo et al. 2000 and references therein) suggest that the manner in which 
the mass stream approaches the disk prevents a hot spot from forming. Instead, a portion 
of the stream makes a full revolution around the disk and then interacts with the original 
stream. The process of this interaction allows the material that has made a full revolution 
around the secondary star to become part of mass stream again; Bisikalo et al. (2000) do not 
consider it to be part of the disk. Since this material's position angle is the same as the disk 
in the system, polarimetry cannot distinguish between the two possibilities. Therefore, we 
use the term 'hot spot' to refer to the region where the mass stream interacts with material 
already encircling the secondary star and assume that any material which has completed a 
revolution around the secondary star is part of the disk. 

Even if a true "hot spot" is not created by the mass stream-disk interaction, the region 
where the stream and disk meet could potentially decrease observed polarization from the 
disk edge by disrupting the otherwise smooth structure of the disk edge and adding unpo- 
larized light at phases when it is visible. On the other hand, the mass stream, which is 
elongated in the same direction as the disk, should produce a polarization position angle 
very similar to that of the disk. Therefore, the presence of the mass stream should not lead 
to a decrease in the observed polarization. The effects of a hot spot would be detectable in 
the polarization light curves in the BVRI bands because the disk is the primary scattering 
region for visible light in the /3 Lyr system. But, if it is not significantly brighter than the 
surrounding disk, the hot spot would not be visible in the total light curves. Therefore, 
we interpret the % Q p minimum associated with the randomization of the polarized position 
angle prior to sec ondary eclipse as the first direct evid ence for the proposed hot spot on the 



f3 Lyr disk edge (ILubow fe Shul 1 19751 ; lHarmaned 120021 ) . 



We expect the hot spot to create an unstructured region on the disk edge where the 
polarization vectors of the scattered light are randomized in position angle. When this part 
of the disk is visible, the hot spot should cause a decreased polarization signal and a rotation 
in position angle, both of which occur in our polarization curves (Figures HJH]). As long 
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as the hot spot does not lie on the line connecting the centers of mass of the two stars, 
and its brightness in the visible continuum is similar to that of the undisturbed disk, its 
effect should result in a minimum in polarization that does not correspond to a minimum 
in flux. Hydrodynamical modeling of the f3 Lyr system indicates that the mass stream, and 
therefore its associated hot spot, should lead the los er in the sense of rotation of the system 



( jLubow fc Shulll975l ; iNazarenko fc Glazunoval 120061 ). In the (3 Lyr polarization curves, the 



polarization minimum and concurrent position angle variation occur just before secondary 
eclipse, suggesting that the hot spot begins its transit of the disk before the loser does. In 
fact, the larger dispersion of points in the B and V bands at the first quadrature phase 
when compared to the second quadrature phase suggests that the hot spot is already in view 
by phase 0.25. In this picture, the minimum polarization occurs at the phases where the 
disk area disrupted by the hot spot and eclipsed by the loser is maximized. We sketch this 
proposed interpretation in Figure H 

We note that in all bands, there are fewer data points after secondary eclipse than before, 
which may skew the %Q P Fourier fit near secondary minimum. However, we have several 
reasons to believe the eclipse offset is not an artifact of the fit. The effect is apparent in all 
filters, some of which have a much lower point density in phase then the V band. (However, 
the V band displays the smallest difference between the %Q P near secondary minimum and 
phase 0.5.) The position angle rotation does not heavily depend on the number of points, 
is apparent in all bands, and has a larger effect at pre-secondary eclipse phases than post- 
secondary eclipse phases. Finally, the uncertainties on the phases at which the minima occur 
are small (see Table [T2"j) compared to the difference between phase 0.5 and the polarization 
secondary minimum. Future work will include filling in the post-secondary data gap with 
new HPOL observations to improve the Fourier fits and quantify the %Q P near secondary 
minimum offsets more reliably. 

In the subsections below, we outline three different estimates of the size of the hot spot, 
assuming it has the same height as the edge of the disk. We use the following values for 
system parameters: a loser radius of Rl = 15R Q , a disk 
separation of Rs = 58i? , and a disk height of H D = 16J 



Linnell 


2000; 


Harmanec 


uiiiai 

2002) 



4.1. Hot Spot Size Estimate: %Q P Method 

We can use the offset in secondary eclipse to estimate a maximum size for the hot spot. 
Assuming circular orbits, we have the scenario depicted in Figure |HJ Knowing that phase 
0.5 occurs at an angle of 180° on the circle depicting the loser's orbit, we can use a simple 
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ratio to find the angle 9, 

^ = ^ a) 

180° 180° -9 w 
where P is the phase for which secondary eclipse occurs in polarized light (Table fT2l) and 
180° — 9 is the angle from zero at which phase P occurs. If we know 9, we can also find the 
length of line x, 

x = R s sm(9) (2) 

where Rs is the radius between the center of the disk and the center of the loser. With the 
length of line x we can estimate the projected size of the hot spot, HSq (hatched region in 
the Observer's View in Figured]), with the following equation, 

r 6Oi? for x > \D D , 

HSq = I (x-R L ) + \D D for \D D > x > R L , (3) 
{ \D D - (R L -x) for R L > x > 

where Dp is the diameter of the accretion disk and Rl is the radius of the loser. 

Using the above formulae we can calculate the maximum projected hot spot size for the 
BVRI bands. Table [121 lists the results. The maximum hot spot size ranges from 22 R Q to 
33 R Q . Since we assume the hot spot has the same height as the disk, these values represent 
'widths' along the projected face of the disk. We do not calculate formal error bars on these 
estimates because the estimates vary so widely. 



4.2. Hot Spot Size Estimate: Position Angle Method 

We also used the variations in position angle to estimate a maximum size for the hot 
spot. First, we calculated the size of the disk in phase. To do this we solved Equation 2 
for 9 when x = \D D . The phase for the left side of the disk as depicted in Figure [8] is then 
given by Equation 1. We calculate this phase to be 0.413. Similarly for the right side of the 
disk we calculate a phase of 0.587. The resulting size of the disk in phase is the difference of 
these phases, or 0.174. 

We then estimated the size of the hot spot in phase by finding how long the random- 
ization of position angle lasts. We assumed any points near secondary eclipse that deviated 
significantly from the average position angle were due to the hot spot. We did not use a 
formal calculation to find these points; we chose the smallest and largest phases for the ran- 
domized position angles by eye. We then took the difference in phase between the deviant 
observations with the smallest and largest phases to calculate the size of the hot spot in 
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phase, 6hs (Table [12]) . In this case the maximum hot spot size, HSpa, is given by the ratio 

(4) 



6Oii! HSp A 



0.174 9 HS 

Table [T2] lists the maximum projected size of the hot spot across the edge of the disk for 
each band. Our hot spot size estimates found with this method range from 26 R Q to 58 R Q . 



4.3. Hot Spot Size Estimate: Simple Model 

We used a simple model for a third estimate of the size of the hot spot. For this model 
we assume the polarization of the disk is uniform across the disk edge. We first calculated 
a baseline qfoc-, the polarized flux due to the disk's self-illumination, by taking the error- 
weighted mean % Q P multiplied by the normalized Fourier fit flux of the observations between 
phases 0.7 and 1.2 for each band (HWN). This assumes all the polarized flux at these phases 
is due to light originating within the disk rather than from the loser, a reasonable assumption 
given the results of HWN. We also define qfmin, the minimum polarized flux near secondary 
eclipse due to the primary star's eclipse and hot spot's transit of the disk, to be the error- 
weighted mean %Q P multiplied by the normalized Fourier fit flux for observations between 
phases 0.4 and 0.55. If we subtract from qfoc the amount of polarized flux blocked by the 
primary star and disrupted by the transit of the hot spot, the result should be qfmin, the 
polarized flux observed at secondary eclipse. 

The amount of polarized flux lost due to the primary eclipsing the disk and the hot spot 
transiting the disk is given by 

n if DC , qfDC , , 
qjDC Z A ecl A H S = qfmin (5) 

Ad a d 

where A ec i is the area of the disk eclipsed by the primary star (shaded region in Figure 
HH), Ad is the observed area of the edge of the disk, and Ahs = HdHSsm is the area of 
the disk edge disrupted by the hot spot. The fraction qfoc I Ad = qfocl^DUo gives the 
polarized flux per unit area from the disk. The second term is the polarized flux eclipsed by 
the primary star and the third term is the polarized flux subtracted by the hot spot. We 
note that Equation 5 assumes an inclination angle of i = 90° and that each unit area of the 
disk contributes to the polarized flux equally. In reality, the relative contributions of each 
portion of a disk with an inclination angle of % = 90° are not equal due to limb darkening. 
Taking limb darkening into account would complicate our estimates; a hot spot near a limb 
darkened edge of the disk would need to be larger to account for the same amount polarized 
flux loss that would be lost by a hot spot closer to the center of the disk. Table [12] gives 
qfoc, qfmin, and the resulting hot spot size estimate for the BVRI bands using Equation 5. 
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Assuming an inclination angle of i = 86° ( ILinnell et al.lll998l ; iLinnelll 120001 ) instead of 
90° changes the equation we use to calculate HSsm- In this case, the area of the disk that 
we see is larger than in the edge-on case; projection effects allow us to see a small portion 
of the back side of the disk. The visible portion of the interior of the disk is polarized 
differently than the disk edge; the polarized flux from the interior should cancel with some 
of the polarized flux from the disk edge. Also, the area eclipsed by the primary star is larger 
and the center of the star is no longer aligned with the plane that cuts the disk into equal 
bottom and top halves. 

In order to calculate an estimate for i = 86°, we make the assumption that the area of 
the disk we see is rectangular. This makes our calculations easier since the projected height 
of the disk does not change across the disk edge, but has the effect of making our estimate 
smaller; if the edges of the projected disk have a height of H D and the center of the projected 
disk has a larger height due to the projection, the total area of the disk is smaller than a 
rectangle whose height is the projected height (see Figure [9]). This assumption allows us to 
use Equation 5 with a slight adjustment: 

r if DC , if DC a r s a \ 

QJDC 77 — 77 A ed — — — — AhS — qjmin (Oj 

where Hpp, is the projected disk height. A simple calculation reveals that Hpp, = 20R Q . 
The area of the hot spot, A HS = H d HSsm, remains unchanged because the hot spot is 
only on the front portion of the disk. Therefore, it does not take up the full projected disk 
height. We note that we do not have to account for the cancellation of polarized light due 
to the contribution from the interior of the disk. We are using our observations to estimate 
if dc an d qfmin, and these numbers should therefore already incorporate this effect, if it is 
present. However, we still make the assumption that each portion of the disk contributes 
to the polarized flux equally. Besides the complication due to limb darkening mentioned 
previously, this enlarges our size estimate because more area is needed to cancel out the 
same amount of polarized flux. Table [12] gives the resulting hot spot size estimate for the 
BVRI bands when % = 86° using Equation 6. 

In our estimates, the size of the hot spot is smaller when the inclination angle is 86° 
compared to 90° for the following reason. The area of the disk is larger in the i = 86° method 
compared to the i = 90° method by a factor of approximately 1.3. This causes the amount 
of polarized flux per unit area to decrease by a factor of 0.8. However, the area eclipsed by 
the primary star increases by more than a factor of two. Therefore, the amount of polarized 
light lost due to the hot spot is smaller when % = 86° than when % = 90°. 
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4.4. Comparison and Review of Hot Spot Size Estimates 

Comparing all three methods, we find a wide range of sizes for the hot spot. The 
smallest size estimate is 2 R Q (I band) while the largest is 58 R (V band). The R and 
/ band estimates are the most likely to change with additional data because their current 
phase coverage is not as good as the B and V bands. The B band estimates have the closest 
agreement between the three methods; they range from 22 R to 33 R , while the V band 
estimates have the largest range, 9 R to 58 R . The overlap region for size ranges in all 
the bands is 22 R to 33 R . 

The large size of our estimates lends support to the possibility that we are actually 
detecting the portion of the mass stream which has not interacted with the disk and not a 
hot spot, similar to the findings of Bisikalo et al. (2000). The largest hot spot size estimate, 
58 R , is similar in size to the diameter of the disk, 60 R , and the same as the binary 
separation, 58 R . However, this scenario would not likely produce the phenomena seen in 
Figures [T] through H] because light scattering from the mass stream would tend to have a 
position angle similar to that of the disk; therefore, we prefer a large hot spot interpretation. 

The 22 R to 33 R range is likely an upper limit for the size of the hot spot for several 
reasons. The position angle method (see Section 4.2) relies on using the randomization 
of position angles around secondary eclipse. We determined the length of time that the 
randomization in position angle lasts by using data from multiple orbits of the system. If 
the spot varies, either in size or location, on time scale similar to the orbital period, this 
variation would cause our estimate to be larger than the actual size of the hot spot. Also, we 
made the assumption that the hot spot contributes only unpolarized light to the observations. 
It is possible that the hot spot contributes light polarized at a different position angle than 
light polarized in the disk. This would cause a cancellation effect to occur; light polarized in 
the hot spot would cancel some of the light polarized in the disk. In this case, our estimate 
would again be larger than the area of the hot spot. 

If the hot spot is indeed larger than 30 R , a portion of it may already be visible to 
the observer during primary eclipse. If this were the case, we would expect the polarization 
maxima in the affected bands to shift from phase 0.0 to an earlier phase. This line of 
reasoning suggests that the primary eclipse effects seen in Figures [1] through HI if real, may 
also be due to the hot spot (see Section 3.1). However, such a scenario does not explain 
why the / band polarization maximum occurs just after primary eclipse. The R band and 
three-term V band Fourier fits do not show an offset at primary eclipse, suggesting that the 
hot spot is not visible at this phase (see Section 3.1). We propose that the primary eclipse 
offsets are likely a result of the Fourier fits in some bands being less well determined around 
primary eclipse due to incomplete phase coverage. 
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The position angle scatter around primary eclipse in Figures [T] through H] could also 
be interpreted as evidence that the hot spot is visible at this phase. As the system moves 
from primary eclipse to secondary eclipse, the amount of the hot spot visible to the observer 
should only increase if the time scale for changes in the hot spot is large compared to the 
orbital period. This suggests that if the scatter in the position angle of the 1997 August 
25 observations (squares between phases 0.000 and 0.015 in Figures [1] through @| is due to 
the hot spot, then a similar scatter should also exist in observations where the hot spot 
was fully visible during the same orbit. In particular, the 1997 August 26 (phase 0.079) 
and 1997 August 30 (phase 0.385) observations should show this effect. However, these two 
observations appear to have a position angle more similar to the system average than the 
1997 August 25 observations. If that the mass transfer rate varies on time scales shorter 
than an orbital period, then the hot spot size may have changed over this five-day period to 
change the amount of position angle scatter. This scenario would explain the position angle 
scatter at primary eclipse for a single cycle, but it would not address the fact that no other 
orbits of the system were observed to have a position angle scatter near primary eclipse. 
The combination of the arguments for and against the visibility of the hot spot at primary 
eclipse does not clearly determine whether its effects are discernible near phase 0.0. Multiple 
spectropolarimetric observations from the same orbit near primary eclipse will provide more 
insight into the cause of this phenomena and access its repeatability. In particular, the 
evolution of the color index of the polarized flux through primary eclipse compared to the 
evolution of the color index of the total light flux may shed new light on this situation. 

Harmanec et al (1996) derived the location of the 'roots' of the bipolar outflows, where 
the bipolar outflows originate within the disk, using the Ha abs orption cor e and the Ha 



emission wings seen in many epochs of spectra (see Figure 1 in lHarmaned (120021 ) for an 
artist's view of the location of the outflows). The location of the roots (marked by a filled 
square in Figure [H] for the Ha absorption and by a filled star for Ha emission) is in the 
same quadrant of the disk in which we interpret the hot spot being located. However, 
the location of the roots is in the interior of the disk while we suggest the hot spot is a 
disruption in the structure of the disk edge. Certainly the bipolar outflows and the hot spot 
are related; both components are the result of the system's high mass transfer rate. How 
far into the interior of the disk the hot spot reaches is unknown. Addition ally, the disk i s 



made up of two components: a dense inner disk and an outer less dense disk ( )Skulskiilll992l ). 
What constitutes the 'disk edge' where the hot spot disruption occurs is unclear, although 
we have assumed it is the outermost edge of the disk (whose diameter is 60 R Q ) in our size 
estimates. Because the scattering region for the Ha line is thought to be the bipolar outflows 
(HNF), future high precision Ha line polarization measurements may be able to link the two 
structures. If the hot spot location is consistent with the roots of the outflows, we expect 
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the Ha line polarization to show a secondary minimum offset, similar to those characterizing 
the broadband curves (Figures [1] through HJ), while maintaining a position angle consistent 
with the outflows. Radial velocity curves of the Ha line's polarized flux may also provide 
valuable insight into the relationship between the hot spot and bipolar outflows. 



5. SUMMARY 

We have presented a large new data set of polarimetric observations of [3 Lyr in the 
BVRI bands and the first Fourier fits to the polarimetric variations in these bands and 
the He I A5876 emission line. We have interpreted the minimum in the BVRI projected 
polarization prior to secondary eclipse and the associated position angle rotations as the 
first direct evidence for a hot spot on the edge of the accretion disk in the (3 Lyr system. 
Using the phases of polarization minimum, the scatter of the position angle and a simple 
model, we have estimated the maximum size of the hot spot to be between 22 and 33 R Q 
across the face of the disk. More extensive polarimetric modeling of (3 Lyr is needed in order 
to fully understand these results. Insights into the importance of the effects at primary 
eclipse and more accurate estimates of the hot spot size could be derived from such models. 

We expect the hot spot may also be detectable in X-rays. Both ROSAT HRI (Berghofer 
& Schmitt 1994) and Suzaku (Ignace et al. 2008) have detected strong and variable hard 
X-ray emission from (3 Lyr. However, neither set of observations has provided information 
on the origin of the X-ray emission or observed the system at phases at which we see the hot 
spot effects. An X-ray light curve with more complete phase coverage will help locate the 
source of the X-ray emission. 

The large uncertainties and scatter in Figure |6]make it difficult to pinpoint a location for 
the origin of the jets with confidence. Future, higher-precision line polarization measurements 
will provide much needed insights and determine their source. 

Advancements in technology will soon allow fo r the combination of long-baseline optical 



interferometry with polarimetry (jElias et al.l 120081 ). We expect such a technological devel- 



opment will provide new and exciting geometrical insights into the (3 Lyr system and others 
like it. 
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Table 1. Date and Phase Information for Midpoints of the HPOL f3 Lyrae and f3 Lyrae B 

Observations 



Date HJD - 2,400,000 Phase a 

j3 Lyr Reticon (through HNF; recalibrated): 



1992 SeD 30 


48895 71 


0.661 


1992 Oct 6 


48901.67 


0.121 


1992 Oct 13 


48908.60 


0.657 


1992 Oct 27 


48922.58 


0.737 


1992 Dec 28 


48984.52 


0.525 


1993 Jul 26 


49194.60 


0.762 


1994 Jun 3 


49506.78 


0.890 


1994 Jun 30 


49533.59 


0.962 


1994 Jul 22 


49555.66 


0.668 


1994 Jul 29 


49562.62 


0.206 


1994 Jul 31 


49564.59 


0.358 


1994 Sep 8 


49603.53 


0.368 


1994 Sep 20 


49615.49 


0.292 


1994 Nov 7 


49663.46 


0.000 


P Lyr CCD (through HNF, recalibrated): 


1995 Mar 14 b 


49790.89 


0.848 


1995 May 5 


49842.82 


0.862 


1995 May 26 


49863.82 


0.485 


1995 May 27 c 


49864.78 


0.559 


1995 May 30 


49867.84 


0.796 


1995 Jun 4 


49872.86 


0.184 


1995 Jul 3 


49901.66 


0.410 


1995 Jul 10 


49908.72 


0.955 


1995 Jul 12 


49910.82 


0.117 


1995 Jul 18 


49916.70 


0.572 


1995 Jul 24 


49922.79 


0.043 


1995 Aug 6 


49935.68 


0.039 


1995 Aug 14 c 


49943.63 


0.653 


1995 Aug 18 


49947.72 


0.969 


1995 Sep 10 


49970.74 


0.749 



Table 1 — Continued 



Date HJD - 2,400,000 Phase a 



p Lyr CCD (after HNF, final calibration): 



1996 Aug 21 


50316. 


.65 


0.483 


1997 May 17 d 


50585. 


.78 


0.283 


1997 May 26 d 


50594. 


.80 


0.980 


1997 Jul 5 


50634. 


.82 


0.073 


1997 Jul 7 


50636. 


.72 


0.219 


1997 Jul 7 


50636. 


,80 


0.226 


1997 Jul 10 


50639. 


.67 


0.447 


1997 Jul 10 


50639. 


.77 


0.455 


1997 Jul 10 


50639. 


.87 


0.463 


1997 Jul 11 


50640. 


.68 


0.525 


1997 Jul 11 


50640. 


.79 


0.534 


1997 Jul 15 


50644. 


.66 


0.833 


1997 Jul 15 


50644. 


.77 


0.842 


1997 Jul 15 


50644. 


.86 


0.848 


1997 Jul 16 


50645. 


.68 


0.912 


1997 Jul 16 


50645. 


.78 


0.920 


1997 Jul 18 


50647. 


71 


0.069 


1997 Aug 1 


50661. 


.64 


0.145 


1997 Aug 25 


50685. 


.64 


0.000 


1997 Aug 25 


50685. 


.74 


0.008 


1997 Aug 25 


50685. 


.83 


0.015 


1997 Aug 26 


50686. 


.66 


0.079 


1997 Aug 30 


50690. 


.62 


0.385 


1997 Sep 7 


50698. 


63 


0.004 


1997 Sep 7 


50698. 


71 


0.010 


1997 Sep 12 


50703. 


.77 


0.401 


1997 Sep 21 


50712. 


.64 


0.087 


1997 Sep 25 


50716. 


.73 


0.403 


1997 Oct 3 


50724. 


60 


0.011 


1997 Oct 3 


50724. 


.70 


0.019 
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Table 1 — Continued 



Date HJD - 2,400,000 Phase a 



1 QQ7 Oct 4 


50725.57 


0.086 


1 QQ7 OH- 98 


50749.54 


0.939 


1997 Nov 17 


50769.56 


0.486 


1997 Dec 15 


50797.50 


0.645 


1998 Apr 19 


50922.86 


0.333 


1998 Apr 24 


50927.90 


0.723 


1998 Jun 23 


50987.82 


0.354 


1998 Jul 31 


51025.78 


0.287 


1998 Aug 31 


51056.62 


0.671 


1998 Sep 8 


51064.76 


0.300 


p Lyr B CCD (throu 


gh HNF, recalibrated): 


1995 May 21 


49858.85 




1996 Jul 3 


50267.75 




f3 Lyr B CCD (after HNF, final calibration): 


1998 Dec 8 C 


51155.54 




1998 Dec ll c 


51158.53 




1998 Dec 12 d 


51159.52 




1999 Aug 30 


51420.63 




1999 Sep 13 


51434.64 




1999 Oct 3 d 


51454.59 




1999 Nov l d 


51483.60 




1999 Nov 2 


51484.57 




1999 Nov 15 d 


51497.55 





a Phases were calculated using the ephemeris 
in Harmanec & Scholz (1993). 

b Table 1 of HNF incorrectly lists this date as 
1994 March 14. 

c These observations used only the red grat- 
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ing; see Section 2.1. 

d These observations used only the blue grat- 
ing; see Section 2.1. 
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Table 2. Date and Phase Information for Midpoints of the FCO (5 Lyrae Observations 



Date HJD - 2,400,000 Phase 



B Band: 








1989 Apr 12 


47628.88 


0. 


.753 


1989 Apr 20 


47636.85 





,352 


1989 Apr 23 


47639.85 





.583 


1989 May 29 


47675.80 





.366 


1989 Jun 5 


47682.78 


0. 


.907 


1989 Jun 12 


47689.80 


0. 


.448 


1989 Jun 19 


47696.78 





,989 


1989 Jun 30 


47707.76 


0. 


.840 


1989 Jul 2 


47709.79 


0. 


,994 


1989 Jul 25 


47732.74 


0. 


,772 


1989 Jul 29 


47736.79 


0. 


,081 


1989 Aug 31 


47769.68 





,632 


1989 Sep 3 


47772.58 





,844 


1989 Sep 4 


47773.57 


0. 


.922 


1989 Sep 5 


47774.57 


0. 


.999 


1989 Sep 9 


47778.57 


0. 


.308 


1989 Oct 5 


47804.45 





.318 


1989 Oct 8 


47807.60 


0. 


,550 


1992 Jun 4 


48777.76 





.543 


V Band: 








1987 Aug 11 


47018.61 





.563 


1987 Aug 13 


47020.57 


0. 


.718 


1987 Aug 15 


47022.59 


0. 


.872 


1987 Aug 16 


47023.58 





.949 


1987 Aug 18 


47025.54 


0. 


.104 


1987 Aug 18 


47025.56 





.104 


1987 Aug 20 


47027.59 


0. 


.259 


1987 Aug 21 


47028.58 





.336 


1987 Aug 22 


47029.56 


0. 


,413 


1987 Aug 24 


47031.58 


0. 


.568 
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Table 2 — Continued 



Date HJD - 2,400,000 Phase a 



1987 


Sep 15 


47053 


.57 





.268 


1987 


Sep 26 


47064 


.55 





.119 


1987 


Sep 27 


47065. 


.54 





.196 


1987 


Sep 28 


47066 


.55 





.273 


1987 


Oct 5 


47073. 


.55 





.814 


1987 


Oct 6 


47074. 


.53 





.892 


1988 


Jun 13 


47325. 


.80 





.313 


1988 


Jun 24 


47336 


.80 





.163 


1988 


Jun 28 


47340. 


.61 





.453 


1988 


Jul 1 


47343 


.62 





.685 


1988 


Jul 16 


47358. 


.66 





.863 


1988 


Aug 19 


47392. 


.69 





.491 


1988 


Aug 22 


47395. 


.68 





.723 


1988 


Aug 23 


47396. 


.59 





.781 


1988 


Sep 1 


47405. 


.62 





.477 


1988 


Sep 2 


47406. 


.58 


0. 


.554 


1988 


Sep 3 


47407. 


.59 





.631 


1988 


Sep 7 


47411. 


.56 


0. 


.941 


1988 


Sep 8 


47412. 


.55 





.018 


1988 


Sep 11 


47415. 


.55 





.250 


1988 


Sep 15 


47419. 


.55 





.559 


1988 


Sep 16 


47420. 


.55 





.636 


1988 


Sep 27 


47431 


.54 





.487 


1988 


Sep 28 


47432. 


.50 





.564 


1988 


Sep 29 


47433. 


.53 





.641 


1988 


Oct 1 


47435. 


.57 


0. 


.796 


1988 


Oct 2 


47436 


.54 





.873 


1988 


Oct 6 


47440. 


.56 


0. 


.182 


1988 


Oct 10 


47444. 


.52 





.491 


1988 


Oct 14 


47448 


.52 





.801 


1988 


Oct 15 


47449. 


.57 





.878 
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Table 2 — Continued 



Date HJD - 2,400,000 Phase a 



1988 


Oct 26 


47460. 


.47 





.728 


1988 


Oct 27 


47461. 


.57 





.805 


1988 


Oct 28 


47462. 


53 


0. 


.883 


1988 


Oct 30 


47464. 


.50 





.037 


1988 


Oct 31 


47465. 


.49 





,115 


1988 


Nov 4 


47469. 


.50 





.424 


1988 


Nov 10 


47475. 


53 





,888 


1988 


Nov 14 


47479. 


.50 





.197 


1988 


Nov 15 


47480. 


.47 





.274 


1988 


Nov 18 


47483. 


.49 


0. 


.506 


1989 


Apr 12 


47628. 


.85 


0. 


.733 


1989 


Apr 20 


47636. 


.83 


0. 


.352 


1989 


Apr 23 


47639. 


.83 





.583 


1989 


May 29 


47675. 


.78 


0. 


.366 


1989 


Jun 5 


47682. 


.80 





.907 


1989 


Jun 12 


47689, 


.82 


0. 


,448 


1989 


Jun 19 


47696. 


.80 


0. 


.989 


1989 


Jun 30 


47707. 


.78 


0. 


.840 


1989 


Jul 2 


47709. 


.81 


0. 


.994 


1989 


Jul 25 


47732. 


.76 





.772 


1989 


Jul 29 


47736. 


.82 





.081 


1989 


Aug 31 


47769. 


.70 





,632 


1989 


Sep 3 


47772. 


60 





,844 


1989 


Sep 4 


47773. 


.61 





.922 


1989 


Sep 5 


47774. 


60 





,999 


1989 


Sep 9 


47778. 


.60 


0. 


.308 


1989 


Oct 2 


47801. 


.68 





.105 


1989 


Oct 23 


47822. 


.59 





.709 


1992 


Jun 13 


48786. 


.72 


0. 


.239 


1992 


Jun 18 


48791. 


.65 





.625 


1992 


Jun 21 


48794. 


.62 





.838 
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Table 2 — Continued 



Date HJD - 2,400,000 Phase a 



1992 Jun 28 


48801. 


73 


0.398 


1992 Jun 29 


48802. 


.67 


0.476 


1992 Jul 6 


48809. 


.68 


0.017 


1992 Jul 8 


48811. 


.59 


0.152 


1992 Jul 19 


48822. 


60 


0.002 


1992 Jul 29 


48832. 


.57 


0.775 


1992 Jul 30 


48833, 


.58 


0.852 


1992 Aug 2 


48836. 


.61 


0.084 


1992 Aug 3 


48837. 


60 


0.161 


1992 Aug 5 


48839. 


60 


0.316 


1992 Aug 8 


48842. 


.56 


0.548 


1992 Aug 11 


48845. 


.56 


0.780 


1992 Aug 19 


48853. 


.54 


0.398 


1992 Aug 30 


48864. 


.56 


0.248 


1992 Sep 1 


48866. 


.60 


0.403 


1992 Sep 13 


48878, 


.55 


0.330 


R Band: 








1989 Apr 12 


47628. 


90 


0.753 


1989 Apr 20 


47636. 


.81 


0.352 


1989 Apr 23 


47639. 


.79 


0.583 


1989 May 29 


47675. 


.76 


0.366 


1989 Jun 5 


47682. 


.83 


0.907 


1989 Jun 12 


47689. 


.84 


0.448 


1989 Jun 19 


47696. 


.82 


0.989 


1989 Jun 30 


47707. 


.80 


0.840 


1989 Jul 2 


47709. 


.83 


0.994 


1989 Jul 25 


47732. 


.78 


0.772 


1989 Jul 29 


47736. 


.84 


0.081 


1989 Aug 31 


47769. 


.72 


0.632 


1989 Sep 3 


47772. 


.62 


0.844 


1989 Sep 5 


47774. 


.62 


0.999 
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Table 2 — Continued 



Date 


HJD - 2,400,000 


Phase a 


1989 Sep 9 


47778.62 


0.308 


1989 Oct 2 


47801.65 


0.105 


1989 Oct 10 


47809.61 


0.704 



a Phases were calculated using the 
ephemeris in Harmanec & Scholz (1993). 



Table 3. FCO B Filter Data 



Phase 


HJD - 2,400,000 


%Q 


%U 


ISP sub %Q 


ISP sub % U 


°i,Q a (%) 


^,c/ a (%) 


0.081 


47736.79 


0.395 


-0.432 


0.164 


-0.067 


0.007 


0.007 


0.308 


47778.57 


0.579 


-0.493 


0.348 


-0.128 


0.016 


0.014 


0.318 


47804.45 


0.617 


-0.527 


0.386 


-0.162 


0.010 


0.008 


0.352 


47636.85 


0.405 


-0.438 


0.174 


-0.073 


0.008 


0.008 


0.366 


47675.80 


0.549 


-0.520 


0.318 


-0.155 


0.009 


0.008 


0.448 


47689.80 


0.257 


-0.261 


0.026 


0.104 


0.007 


0.007 


0.543 


48777.76 


0.148 


-0.350 


-0.083 


0.015 


0.009 


0.009 


0.550 


47807.60 


0.460 


-0.698 


0.229 


-0.333 


0.016 


0.023 


0.583 


47639.85 


0.421 


-0.527 


0.190 


-0.162 


0.008 


0.010 


0.632 


47769.68 


0.497 


-0.517 


0.266 


-0.152 


0.010 


0.010 


0.753 


47628.88 


0.538 


-0.479 


0.307 


-0.114 


0.014 


0.012 


0.772 


47732.74 


0.486 


-0.494 


0.255 


-0.129 


0.011 


0.011 


0.840 


47707.76 


0.416 


-0.469 


0.185 


-0.104 


0.007 


0.008 


0.844 


47772.58 


0.406 


-0.476 


0.175 


-0.111 


0.009 


0.011 


0.907 


47682.78 


0.531 


-0.488 


0.300 


-0.123 


0.012 


0.011 


0.922 


47773.57 


0.470 


-0.456 


0.239 


-0.091 


0.010 


0.010 


0.989 


47696.78 


0.602 


-0.566 


0.371 


-0.201 


0.010 


0.010 


0.994 


47709.79 


0.635 


-0.558 


0.404 


-0.193 


0.018 


0.016 


0.999 


47774.57 


0.553 


-0.645 


0.322 


-0.280 


0.009 


0.010 



a The CTi Q and j,^ columns represent internal errors. 
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Table 4. FCO V Filter Data 



Phase HJD - 2,400,000 %Q %U ISP sub %Q ISP sub %t/ o"i iQ a (%) a ijC/ a (%) 



0.002 


48822.60 


0.353 


-0.393 


0.134 


-0.028 


0.053 


0.053 


0.017 


48809.68 


0.405 


-0.471 


0.186 


-0.106 


0.013 


0.013 


0.018 


47412.55 


0.508 


-0.433 


0.289 


-0.068 


0.010 


0.009 


0.037 


47464.50 


0.687 


-0.704 


0.468 


-0.339 


0.097 


0.099 


0.081 


47736.82 


0.390 


-0.422 


0.171 


-0.057 


0.010 


0.011 


0.084 


48836.61 


0.390 


-0.423 


0.171 


-0.058 


0.010 


0.010 


0.104 


47025.54 


0.287 


-0.412 


0.068 


-0.047 


0.014 


0.019 


0.104 


47025.56 


0.292 


-0.437 


0.073 


-0.072 


0.008 


0.012 


0.105 


47801.68 


0.548 


-0.458 


0.329 


-0.093 


0.010 


0.009 


0.115 


47465.49 


0.310 


-0.418 


0.091 


-0.053 


0.006 


0.008 


0.119 


47064.55 


0.335 


-0.454 


0.116 


-0.089 


0.008 


0.011 


0.152 


48811.59 


0.440 


-0.416 


0.221 


-0.051 


0.012 


0.012 


0.161 


48837.60 


0.442 


-0.438 


0.223 


-0.073 


0.010 


0.010 


0.163 


47336.80 


0.370 


-0.390 


0.151 


-0.025 


0.092 


0.096 


0.182 


47440.56 


0.519 


-0.431 


0.300 


-0.066 


0.012 


0.010 


0.196 


47065.54 


0.397 


-0.457 


0.178 


-0.092 


0.009 


0.010 


0.197 


47479.50 


0.451 


-0.474 


0.232 


-0.109 


0.015 


0.016 


0.239 


48786.72 


0.322 


-0.454 


0.103 


-0.089 


0.015 


0.015 


0.248 


48864.56 


0.500 


-0.510 


0.281 


-0.145 


0.006 


0.006 


0.250 


47415.55 


0.492 


-0.466 


0.273 


-0.101 


0.014 


0.013 


0.259 


47027.59 


0.412 


-0.441 


0.193 


-0.076 


0.009 


0.010 


0.268 


47053.57 


0.420 


-0.433 


0.201 


-0.068 


0.007 


0.007 


0.273 


47066.55 


0.493 


-0.448 


0.274 


-0.083 


0.008 


0.008 


0.274 


47480.47 


0.549 


-0.570 


0.330 


-0.205 


0.011 


0.011 


0.308 


47778.60 


0.560 


-0.461 


0.341 


-0.096 


0.013 


0.011 


0.313 


47325.80 


0.408 


-0.419 


0.189 


-0.054 


0.009 


0.009 


0.316 


48839.60 


0.526 


-0.492 


0.307 


-0.127 


0.008 


0.008 


0.330 


48878.55 


0.580 


-0.475 


0.361 


-0.110 


0.012 


0.012 


0.336 


47028.58 


0.461 


-0.423 


0.242 


-0.058 


0.008 


0.008 


0.352 


47636.83 


0.391 


-0.421 


0.172 


-0.056 


0.008 


0.009 


0.366 


47675.78 


0.556 


-0.495 


0.337 


-0.130 


0.012 


0.011 


0.398 


48801.73 


0.455 


-0.428 


0.236 


-0.063 


0.021 


0.021 


0.398 


48853.54 


0.332 


-0.397 


0.113 


-0.032 


0.012 


0.012 


0.403 


48866.60 


0.492 


-0.401 


0.273 


-0.036 


0.009 


0.009 


0.413 


47029.56 


0.312 


-0.268 


0.093 


0.097 


0.044 


0.039 


0.424 


47469.50 


0.405 


-0.414 


0.186 


-0.049 


0.008 


0.008 


0.448 


47689.82 


0.269 


-0.273 


0.050 


0.092 


0.007 


0.007 


0.453 


47340.61 


0.318 


-0.282 


0.099 


0.083 


0.009 


0.008 


0.476 


48802.67 


0.233 


-0.198 


0.014 


0.167 


0.017 


0.017 


0.477 


47405.62 


0.328 


-0.268 


0.109 


0.097 


0.006 


0.005 


0.487 


47431.54 


0.430 


-0.413 


0.211 


-0.048 


0.008 


0.007 


0.491 


47392.69 


0.281 


-0.349 


0.062 


0.016 


0.009 


0.011 


0.491 


47444.52 


0.434 


-0.432 


0.215 


-0.067 


0.013 


0.010 


0.506 


47483.49 


0.501 


-0.401 


0.282 


-0.036 


0.011 


0.009 


0.548 


48842.56 


0.316 


-0.455 


0.097 


-0.090 


0.019 


0.019 


0.554 


47406.58 


0.354 


-0.345 


0.135 


0.020 


0.010 


0.009 



-33 - 



Table 4 — Continued 



Phase 


HJD - 2,400,000 


%Q 


%U 


ISP sub %Q 


ISP sub % U 


<7i,Q a (%) 


CTi,C7 a (%) 


o ^^q 


4741 Q ^ 


n 341 

U.04:_L 


4fi8 

~U.4:UO 


122 


103 
-U. -LUO 


007 
U.UU 1 


OOQ 

u.uuy 


o ^fi3 

U. OOO 


4701 8 R1 
4: ( UlO.Dl 


n 3^ 

U.OOO 


38Q 

-u.ooy 


1 3fi 


094 

-U.UZ4: 


007 
U.UU i 


O 008 
U.UUo 


U.OU4: 


47439 ^0 

4: I 4:OZ.OU 


n 379 

U.O ( Z 


44^ 

-U.4: 4:0 


U. -LOO 


080 

-u.uou 


01 3 
U . U -L O 


01 
U.UIO 


o ^fi8 

U.OOO 


47031 ^8 

4: I UOl.OO 


n 399 


409 

~U.4:UZ 


1 03 
U. 1UO 


037 
-U.UO 1 


OOfi 

u.uuo 


O 008 
U.UUO 


n ^83 
u. ooo 


47fi3Q 83 


421 


^01 
-U.OUl 


202 


-U. -LOO 


01 3 
U . U -L O 


01^ 
U.UIO 


o 

u.uzo 


487Q1 fi 1 ^ 


0^17 
U.Ol l 


^93 

-u.ozo 


9Q8 

u.zyo 


-U. -LOO 


004 

U.UU4: 


004 

U.UU4: 


n <^3i 


47407 

4: I 4:U < .Oy 


n ^i 4 

U.O-L4: 


438 

~U.4:00 


9Q^ 
u.zyo 


073 
-U.U ( o 


01 1 

U.U-L 1 


OOQ 

u.uuy 


n £39 


477RQ 70 

4: i ( uy. i u 


n 47^ 

U.4: ( O 


4RQ 

-u.4:uy 


9^fi 
u.zoo 


1 04 

-U. -LU4: 


019 
U.U-LZ 


011 
U.U1 -L 


o fi^fi 
u.uou 


47490 ^ 

4: I 4:ZU.OO 


n 398 


947 

"U.Z4: I 


1 OQ 
u. iuy 


118 

U. -L IO 


00Q 

u.uuy 


007 
U.UU / 


U.U4:± 


47433 ^3 


n ^7^ 

U.O I o 


3Q0 

-u.oyu 


U.OOO 


09^ 

-u.uzo 


01 3 
U . U -L O 


010 
U.U1U 


n £8^ 

U.OoO 


4: / 040.DZ 


n ^07 

U. OU I 


438 


988 
U.Zoo 


073 
-U.U I O 


010 
U.U1U 


O OOQ 

u.uuy 


70Q 

u. i uy 


47899 

4: / ozz.oy 


n ^08 

U. OUo 


488 


98Q 

u.zoy 


1 93 
-U. 1ZO 


OOQ 

u.uuy 


O 008 
U.UUo 


U. < io 


47090 ^7 

4: I UZU.O ( 


n ^91 

U. OZ_L 


3Q8 

-u.oyo 


309 

u.ouz 


033 

-u.uoo 


019 
U.U-LZ 


010 
U.U1U 


o 793 

U. f zo 


4730^ R8 

4: / oyo.uo 


n ^^8 


473 

-U.4: I O 


33Q 

u.ooy 


1 08 
-U. -LUO 


01 1 

U.U-L 1 


OOQ 

u.uuy 


o 798 

U. ( ZO 


474R0 47 

4: / 4:0U.4: 1 


478 

U.4: I o 


470 

-U.4: 1 U 


9^Q 

u.zoy 


-U. 1UO 


OOQ 

u.uuy 


O OOQ 

u.uuy 


o 733 

U. f oo 


47(^98 8^ 

4: I OZO.OO 


n ^i 8 

U.OIO 


4^0 

-u.4:oy 


9QQ 
u.zyy 


0Q4 

-u.uy4: 


OOQ 

u.uuy 


008 

u.uuo 


779 


47739 7fi 

4: I I OZ. / O 


480 

U.4:OU 


-0 449 


9fi1 


084 

-U.U04: 


01 1 

U.U-L 1 


010 
U.U1U 


0.775 


48839 ^7 

4:OOOZ. O ( 


0.508 


-0.491 


0.289 


-0.126 


0.011 


0.011 


n 7sn 

U. f ou 


4884^ ^fi 

4:004:0.00 


n 4^8 

U.4:00 


477 

~U.4:( I 


93Q 

u.zoy 


-0 112 


010 
U.U-LU 


010 
U.U1U 


U. ( ol 


473Q£ 

4: ( oyo.oy 


n 483 


41 Q 

-u.4:iy 


9(^4 


0^4 

-U.U04: 


OOIO 
U.U1U 


O OOQ 

u.uuy 


u. * yo 


4: / 4:O0.0 1 


fl 4R9 

U.4:0Z 


417 

~U.4:1 I 


943 

U.Z4:0 


0^9 

-u.uoz 


OOIO 
U.U1U 


O OOQ 

u.uuy 


n 801 


47448 ^9 

4: / 4:4:0. OZ 


499 

U.4:ZZ 


44R 

-U.4: 4:0 


903 
U.ZUO 


081 
-U.Uol 


O 007 
U.UU i 


O 007 
U.UU / 


n 80^ 

U.oUO 


474R1 ^7 

4: I 4:01.0 ( 


fl 491 

U.4:Z_L 


494 

~U.4:Z4: 


909 
U.ZUZ 


O^Q 

-u.uoy 


OOIO 
U.U1U 


OOIO 
U.U1U 


o 81 4 

U.O-L4: 


47073 55 


43^ 

U.4:00 


403 

~U.4:UO 


91 
U.ZIO 


038 

-u.uoo 


OOQ 

u.uuy 


008 

u.uuo 


U.OOO 


487Q4 R9 
4:0 ( y4:.0Z 


n ^03 

U. OUo 


4^fi 

-U.4:00 


984 

U.Z04: 


OQ1 

-u.uy i 


O OOfi 

u.uuo 


O OOfi 

u.uuo 


0.840 


47707.78 


0.411 


-0.410 


0.192 


-0.045 


0.007 


0.007 


0.844 


47772.60 


0.376 


-0.425 


0.157 


-0.060 


0.008 


0.009 


0.852 


48833.58 


0.391 


-0.432 


0.172 


-0.067 


0.011 


0.011 


0.863 


47358.66 


0.435 


-0.357 


0.216 


0.008 


0.009 


0.008 


0.872 


47022.59 


0.451 


-0.640 


0.232 


-0.275 


0.011 


0.010 


0.873 


47436.54 


0.467 


-0.437 


0.248 


-0.072 


0.012 


0.011 


0.878 


47449.57 


0.374 


-0.388 


0.155 


-0.023 


0.009 


0.009 


0.883 


47462.53 


0.389 


-0.397 


0.170 


-0.032 


0.007 


0.008 


0.888 


47475.53 


0.308 


-0.376 


0.089 


-0.011 


0.006 


0.007 


0.892 


47074.53 


0.376 


-0.396 


0.157 


-0.031 


0.010 


0.011 


0.907 


47682.80 


0.468 


-0.442 


0.249 


-0.077 


0.011 


0.011 


0.922 


47773.61 


0.451 


-0.407 


0.232 


-0.042 


0.008 


0.008 


0.941 


47411.56 


0.547 


-0.396 


0.328 


-0.031 


0.013 


0.010 


0.949 


47023.58 


0.540 


-0.490 


0.321 


-0.125 


0.010 


0.009 


0.989 


47696.80 


0.558 


-0.509 


0.339 


-0.144 


0.016 


0.015 


0.994 


47709.81 


0.623 


-0.496 


0.404 


-0.131 


0.016 


0.013 


0.999 


47774.60 


0.514 


-0.603 


0.295 


-0.238 


0.012 


0.014 



a The cr it Q and cr^y columns represent internal errors. 
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Table 5. FCO R Filter Data 



Phase 


HJD - 2,400,000 


%Q 


%U 


ISP sub %Q 


ISP sub % U 


»i,Q B (%) 


^,u a (%) 


0.081 


47736.84 


0.318 


-0.373 


0.131 


-0.077 


0.007 


0.008 


0.105 


47801.65 


0.454 


-0.432 


0.267 


-0.136 


0.011 


0.011 


0.308 


47778.62 


0.503 


-0.431 


0.316 


-0.135 


0.008 


0.007 


0.352 


47636.81 


0.376 


-0.374 


0.189 


-0.078 


0.012 


0.012 


0.366 


47675.76 


0.504 


-0.489 


0.317 


-0.193 


0.018 


0.017 


0.448 


47689.84 


0.261 


-0.267 


0.074 


0.029 


0.006 


0.007 


0.583 


47639.79 


0.377 


-0.450 


0.190 


-0.154 


0.009 


0.011 


0.632 


47769.72 


0.412 


-0.443 


0.225 


-0.147 


0.013 


0.014 


0.704 


47809.61 


0.830 


-0.779 


0.643 


-0.483 


0.017 


0.016 


0.753 


47628.90 


0.452 


-0.417 


0.265 


-0.121 


0.014 


0.013 


0.772 


47732.78 


0.398 


-0.398 


0.211 


-0.102 


0.009 


0.009 


0.840 


47707.80 


0.359 


-0.373 


0.172 


-0.077 


0.009 


0.004 


0.844 


47772.62 


0.330 


-0.372 


0.143 


-0.076 


0.007 


0.007 


0.907 


47682.83 


0.391 


-0.443 


0.204 


-0.147 


0.012 


0.021 


0.989 


47696.82 


0.449 


-0.417 


0.262 


-0.121 


0.011 


0.010 


0.994 


47709.83 


0.524 


-0.438 


0.337 


-0.142 


0.017 


0.015 


0.999 


47774.62 


0.463 


-0.522 


0.276 


-0.226 


0.013 


0.014 



a The <Ti t Q and ct^q columns represent internal errors. 
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Table 6. HPOL Synthetic B Filter Data 



Phase HJD - 2,400,000 %Q %U ISP sub %Q ISP sub %U (%) cr s a (%) 



B Band Rcticon: 












0.000 


49663 


.46 





4323 


-0 


.6592 


0.121 


48901 


.67 


0. 


.3905 


-0. 


.4928 


0.206 


49562. 


.62 





4691 


-0 


.4456 


0.292 


49615 


.49 





.4248 


-0 


.5536 


0.358 


49564 


.59 


0. 


.4874 


-0. 


.4409 


0.368 


49603 


.53 





.4823 


-0 


.4503 


0.525 


48984 


.52 





.2897 


-0 


.4828 


0.657 


48908 


.60 


0. 


.5574 


-0. 


5334 


0.661 


48895 


.71 





.5392 


-0 


.5076 


0.668 


49555 


.66 





.5611 


-0 


.5302 


0.737 


48922 


.58 


0. 


.4881 


-0. 


.5021 


0.762 


49194 


.60 





.4841 


-0 


.5642 


0.890 


49506 


.78 





.4337 


-0 


.4554 


0.962 


49533. 


.59 


0. 


.5317 


-0. 


.6448 


B Band CCD: 












0.000 


50685 


.64 





.3794 


-0 


.8971 


0.004 


50698 


.63 


0. 


.5419 


-0. 


.4592 


0.008 


50685 


.74 





.3222 


-0 


8462 


0.010 


50698 


.71 





.4679 


-0 


.4787 


0.011 


50724. 


.60 


0. 


.4846 


-0. 


.5516 


0.015 


50685 


.83 


0. 


.2608 


-0. 


.8606 


0.019 


50724. 


.70 


0. 


.3952 


-0. 


.5281 


0.039 


49935 


.68 





.3981 


-0 


.6422 


0.043 


49922 


.79 


0. 


.6360 


-0. 


.5666 


0.069 


50647 


.71 


0. 


.4896 


-0. 


.5155 


0.073 


50634 


.82 





.3301 


-0 


.3861 


0.079 


50686. 


.66 


0. 


.3372 


-0. 


.5482 


0.086 


50725. 


.57 


0. 


.2982 


-0. 


.3927 


0.087 


50712 


.64 


0. 


.3415 


-0 


.5013 


0.117 


49910 


.82 





.3875 


-0. 


.4219 


0.145 


50661 


.64 





3608 


-0. 


.4670 


0.184 


49872 


.86 





.2940 


-0 


.4482 


0.219 


50636. 


.72 


0. 


.4525 


-0. 


.5404 


0.226 


50636. 


.80 


0. 


.4000 


-0. 


.5080 


0.283 


50585 


.78 





.4840 


-0 


.6685 


0.287 


51025 


.78 





.3671 


-0 


5343 


0.300 


51064 


.76 





3961 


-0. 


.6073 


0.333 


50922 


.86 





.5179 


-0 


.5637 


0.354 


50987 


.82 





.4517 


-0 


.5322 


0.385 


50690 


.62 


0. 


.4775 


-0. 


.5832 


0.401 


50703. 


.77 


0. 


.3079 


-0. 


.4730 


0.403 


50716 


.73 





.3754 


-0. 


.4673 


0.410 


49901 


.66 





.3475 


-0 


.3311 


0.447 


50639. 


.67 


0. 


.2957 


-0 


3379 


0.455 


50639 


.77 





.3165 


-0 


.3542 



0.2020 


-0.2953 


0.0019 


0.0200 


0.1603 


-0.1289 


0.0030 


0.0200 


0.2388 


-0.0816 


0.0062 


0.0200 


0.1945 


-0.1897 


0.0015 


0.0200 


0.2572 


-0.0770 


0.0023 


0.0200 


0.2521 


-0.0864 


0.0016 


0.0200 


0.0594 


-0.1188 


0.0034 


0.0200 


0.3272 


-0.1695 


0.0030 


0.0200 


0.3090 


-0.1436 


0.0023 


0.0200 


0.3308 


-0.1663 


0.0015 


0.0200 


0.2578 


-0.1382 


0.0019 


0.0200 


0.2538 


-0.2003 


0.0011 


0.0200 


0.2027 


-0.0953 


0.0033 


0.0200 


0.3015 


-0.2810 


0.0019 


0.0200 


0.1492 


-0.5332 


0.0044 


0.0220 


0.3117 


-0.0953 


0.0031 


0.0220 


0.0920 


-0.4823 


0.0046 


0.0220 


0.2377 


-0.1149 


0.0063 


0.0220 


0.2544 


-0.1877 


0.0035 


0.0220 


0.0326 


-0.4982 


0.0082 


0.0220 


0.1650 


-0.1643 


0.0053 


0.0220 


0.1679 


-0.2783 


0.0030 


0.0100 


0.4058 


-0.2028 


0.0034 


0.0100 


0.2594 


-0.1516 


0.0036 


0.0140 


0.0999 


-0.0223 


0.0024 


0.0140 


0.1070 


-0.1843 


0.0029 


0.0220 


0.0680 


-0.0289 


0.0033 


0.0220 


0.1113 


-0.1374 


0.0029 


0.0220 


0.1573 


-0.0581 


0.0023 


0.0100 


0.1307 


-0.1030 


0.0027 


0.0140 


0.0638 


-0.0844 


0.0027 


0.0100 


0.2223 


-0.1765 


0.0028 


0.0140 


0.1698 


-0.1442 


0.0024 


0.0140 


0.2538 


-0.3047 


0.0018 


0.0140 


0.1368 


-0.1704 


0.0028 


0.0150 


0.1659 


-0.2435 


0.0054 


0.0150 


0.2877 


-0.1998 


0.0025 


0.0080 


0.2215 


-0.1683 


0.0025 


0.0150 


0.2473 


-0.2193 


0.0031 


0.0220 


0.0777 


-0.1091 


0.0029 


0.0220 


0.1453 


-0.1027 


0.0057 


0.0220 


0.1174 


0.0327 


0.0024 


0.0100 


0.0655 


0.0261 


0.0025 


0.0140 


0.0863 


0.0097 


0.0028 


0.0140 
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Table 6 — Continued 



Phase 


HJD - 2,400,000 


%Q 


%U 


ISP sub %Q 


ISP sub % U 


<T S a (%) 


<T S * (%) 


0.463 


50639.87 


0.3221 


-0.4072 


0.0920 


-0.0433 


0.0047 


0.0140 


0.483 


50316.65 


0.4778 


-0.4378 


0.2475 


-0.0739 


0.0020 


0.0300 


0.485 


49863.82 


0.3250 


-0.3518 


0.0948 


0.0121 


0.0022 


0.0100 


0.486 


50769.56 


0.3433 


-0.4030 


0.1131 


-0.0391 


0.0041 


0.0220 


0.525 


50640.68 


0.3457 


-0.5121 


0.1155 


-0.1482 


0.0029 


0.0140 


0.534 


50640.79 


0.3062 


-0.5369 


0.0760 


-0.1730 


0.0099 


0.0140 


0.572 


49916.70 


0.4309 


-0.4738 


0.2006 


-0.1099 


0.0020 


0.0100 


0.645 


50797.50 


0.4821 


-0.5407 


0.2519 


-0.1769 


0.0050 


0.0220 


0.671 


51056.62 


0.5691 


-0.5555 


0.3388 


-0.1916 


0.0024 


0.0150 


0.723 


50927.90 


0.5293 


-0.6009 


0.2990 


-0.2369 


0.0051 


0.0080 


0.749 


49970.74 


0.5268 


-0.4740 


0.2965 


-0.1101 


0.0020 


0.0100 


0.796 


49867.84 


0.4758 


-0.5419 


0.2455 


-0.1780 


0.0017 


0.0100 


0.833 


50644.66 


0.3884 


-0.5293 


0.1581 


-0.1654 


0.0024 


0.0140 


0.842 


50644.77 


0.3740 


-0.5620 


0.1438 


-0.1981 


0.0026 


0.0140 


0.848 


49790.89 


0.4561 


-0.5198 


0.2259 


-0.1559 


0.0033 


0.0100 


0.848 


50644.86 


0.3216 


-0.5164 


0.0914 


-0.1526 


0.0025 


0.0140 


0.862 


49842.82 


0.4337 


-0.4914 


0.2035 


-0.1276 


0.0023 


0.0100 


0.912 


50645.68 


0.3993 


-0.4460 


0.1691 


-0.0821 


0.0024 


0.0140 


0.920 


50645.78 


0.4319 


-0.4887 


0.2017 


-0.1249 


0.0041 


0.0140 


0.939 


50749.54 


0.5136 


-0.7251 


0.2834 


-0.3612 


0.0035 


0.0220 


0.955 


49908.72 


0.7422 


-0.5796 


0.5120 


-0.2158 


0.0051 


0.0100 


0.969 


49947.72 


0.6624 


-0.5345 


0.4322 


-0.1707 


0.0041 


0.0100 


0.980 


50594.80 


0.4666 


-0.7196 


0.2365 


-0.3557 


0.0082 


0.0140 



a The <Ji and a s columns represent internal and systematic errors respectively. 
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Table 7. HPOL Synthetic V Filter Data 



Phase HJD - 2,400,000 %Q %U ISP sub %Q ISP sub %U (%) cr s a (%) 



V Band Rcticon: 












0.000 


49663 


.46 





.4485 


-0 


.6204 


0.121 


48901 


.67 


0. 


.3620 


-0. 


.4839 


0.206 


49562 


.62 





.4578 


-0 


.3863 


0.292 


49615 


.49 





.4058 


-0 


.5298 


0.358 


49564 


.59 


0. 


.4499 


-0. 


.4464 


0.368 


49603 


.53 





.4474 


-0 


.4376 


0.525 


48984 


.52 





.2534 


-0 


.4453 


0.657 


48908 


.60 


0. 


.5430 


-0. 


.4910 


0.661 


48895 


.71 


0. 


.4934 


-0 


.4758 


0.668 


49555 


.66 





5306 


-0 


.4902 


0.737 


48922 


.58 


0. 


.4451 


-0. 


.4905 


0.762 


49194 


.60 





.4531 


-0 


.5427 


0.890 


49506 


.78 





.3825 


-0 


.4272 


0.962 


49533. 


.59 


0. 


.5339 


-0. 


.6376 


V Band CCD: 












0.000 


50685 


.64 





3923 


-0 


.7974 


0.004 


50698 


.63 


0. 


.5482 


-0. 


.4414 


0.008 


50685 


.74 





.3562 


-0. 


.7982 


0.010 


50698 


.71 





.5082 


-0 


.4602 


0.011 


50724. 


.60 


0. 


.5266 


-0. 


.5742 


0.015 


50685 


.83 


0. 


.3205 


-0. 


.8314 


0.019 


50724. 


.70 


0. 


.4657 


-0. 


.5516 


0.039 


49935 


.68 





.4228 


-0 


.6091 


0.043 


49922 


.79 


0. 


.6145 


-0. 


.5334 


0.069 


50647 


.71 


0. 


.4805 


-0. 


5109 


0.073 


50634 


.82 





.3706 


-0 


.3841 


0.079 


50686. 


.66 


0. 


.3369 


-0. 


.5427 


0.086 


50725. 


.57 


0. 


.3240 


-0. 


.4117 


0.087 


50712 


.64 


0. 


.3497 


-0 


.4760 


0.117 


49910 


.82 





.3706 


-0. 


.3946 


0.145 


50661 


.64 





3556 


-0. 


.4517 


0.184 


49872 


.86 





3099 


-0 


.4214 


0.219 


50636. 


.72 


0. 


.4250 


-0. 


.5151 


0.226 


50636. 


.80 


0. 


.3969 


-0. 


.5070 


0.283 


50585 


.78 





.4976 


-0 


.6274 


0.287 


51025 


.78 





3863 


-0 


.5250 


0.300 


51064 


.76 





.3968 


-0. 


.5232 


0.333 


50922 


.86 





.5109 


-0 


.5433 


0.354 


50987 


.82 





.4459 


-0 


.4969 


0.385 


50690 


.62 


0. 


.4809 


-0. 


.5536 


0.401 


50703. 


.77 


0. 


.3649 


-0. 


.4747 


0.403 


50716 


.73 





.3778 


-0. 


.4680 


0.410 


49901 


.66 





3568 


-0 


.3198 


0.447 


50639. 


.67 





.3001 


-0 


3333 


0.455 


50639 


.77 





.3256 


-0 


.3258 



0.2292 


-0.2737 


0.0021 


0.0200 


0.1427 


-0.1373 


0.0025 


0.0200 


0.2386 


-0.0398 


0.0056 


0.0200 


0.1865 


-0.1830 


0.0018 


0.0200 


0.2306 


-0.0998 


0.0024 


0.0200 


0.2280 


-0.0909 


0.0020 


0.0200 


0.0339 


-0.0985 


0.0026 


0.0200 


0.3236 


-0.1444 


0.0032 


0.0200 


0.2741 


-0.1291 


0.0021 


0.0200 


0.3113 


-0.1435 


0.0017 


0.0200 


0.2258 


-0.1438 


0.0018 


0.0200 


0.2337 


-0.1960 


0.0013 


0.0200 


0.1590 


-0.0808 


0.0053 


0.0200 


0.3147 


-0.2910 


0.0023 


0.0200 


0.1732 


-0.4508 


0.0029 


0.0100 


0.3288 


-0.0948 


0.0020 


0.0100 


0.1370 


-0.4516 


0.0030 


0.0100 


0.2892 


-0.1139 


0.0037 


0.0100 


0.3075 


-0.2280 


0.0022 


0.0100 


0.1022 


-0.4861 


0.0044 


0.0100 


0.2464 


-0.2053 


0.0029 


0.0100 


0.2038 


-0.2631 


0.0021 


0.0050 


0.3948 


-0.1860 


0.0024 


0.0050 


0.2612 


-0.1644 


0.0024 


0.0070 


0.1516 


-0.0382 


0.0016 


0.0070 


0.1179 


-0.1966 


0.0018 


0.0100 


0.1047 


-0.0652 


0.0021 


0.0100 


0.1304 


-0.1295 


0.0018 


0.0100 


0.1517 


-0.0479 


0.0015 


0.0050 


0.1361 


-0.1047 


0.0019 


0.0070 


0.0903 


-0.0756 


0.0017 


0.0050 


0.2058 


-0.1685 


0.0019 


0.0070 


0.1777 


-0.1604 


0.0016 


0.0070 


0.2777 


-0.2799 


0.0012 


0.0070 


0.1671 


-0.1785 


0.0017 


0.0100 


0.1779 


-0.1771 


0.0030 


0.0100 


0.2916 


-0.1967 


0.0016 


0.0050 


0.2268 


-0.1506 


0.0016 


0.0100 


0.2617 


-0.2070 


0.0019 


0.0100 


0.1452 


-0.1273 


0.0019 


0.0100 


0.1584 


-0.1217 


0.0032 


0.0100 


0.1376 


0.0266 


0.0018 


0.0050 


0.0811 


0.0132 


0.0017 


0.0070 


0.1063 


0.0209 


0.0018 


0.0070 
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Table 7— Continued 



Phase HJD - 2,400,000 %Q %U ISP sub %Q ISP sub %U (%) cr s a (%) 



0.463 


50639.87 


0.3353 


-0.3521 


0.1161 


-0.0056 


0.0031 


0.0070 


0.483 


50316.65 


0.4529 


-0.4384 


0.2338 


-0.0928 


0.0014 


0.0200 


0.485 


49863.82 


0.3189 


-0.3351 


0.1000 


0.0106 


0.0015 


0.0050 


0.486 


50769.56 


0.3571 


-0.4099 


0.1378 


-0.0634 


0.0025 


0.0100 


0.525 


50640.68 


0.3109 


-0.5088 


0.0915 


-0.1620 


0.0019 


0.0070 


0.534 


50640.79 


0.2687 


-0.5048 


0.0508 


-0.1606 


0.0059 


0.0070 


0.559 


49864.78 


0.5740 


-0.4877 


0.3675 


-0.1488 


0.0045 


0.0050 


0.572 


49916.70 


0.4404 


-0.4434 


0.2215 


-0.0974 


0.0014 


0.0050 


0.645 


50797.50 


0.4919 


-0.5245 


0.2736 


-0.1783 


0.0027 


0.0100 


0.653 


49943.63 


0.6036 


-0.4925 


0.3022 


-0.1452 


0.0022 


0.0050 


0.671 


51056.62 


0.5323 


-0.5334 


0.3130 


-0.1869 


0.0015 


0.0100 


0.723 


50927.90 


0.4912 


-0.5611 


0.2714 


-0.2138 


0.0034 


0.0050 


0.749 


49970.74 


0.5054 


-0.4378 


0.2861 


-0.0908 


0.0013 


0.0050 


0.796 


49867.84 


0.4662 


-0.5150 


0.2471 


-0.1688 


0.0012 


0.0050 


0.833 


50644.66 


0.3579 


-0.5014 


0.1386 


-0.1549 


0.0017 


0.0070 


0.842 


50644.77 


0.3457 


-0.5039 


0.1264 


-0.1573 


0.0018 


0.0070 


0.848 


49790.89 


0.4292 


-0.4756 


0.2104 


-0.1292 


0.0020 


0.0050 


0.848 


50644.86 


0.3131 


-0.4951 


0.0939 


-0.1483 


0.0016 


0.0070 


0.862 


49842.82 


0.4222 


-0.4591 


0.2037 


-0.1134 


0.0015 


0.0050 


0.912 


50645.68 


0.3943 


-0.4385 


0.1748 


-0.0918 


0.0016 


0.0070 


0.920 


50645.78 


0.3880 


-0.4313 


0.1690 


-0.0849 


0.0025 


0.0070 


0.939 


50749.54 


0.4909 


-0.6709 


0.2715 


-0.3242 


0.0022 


0.0100 


0.955 


49908.72 


0.6986 


-0.5220 


0.4790 


-0.1747 


0.0035 


0.0050 


0.969 


49947.72 


0.6168 


-0.4929 


0.3989 


-0.1476 


0.0028 


0.0050 


0.980 


50594.80 


0.4644 


-0.5950 


0.2447 


-0.2476 


0.0053 


0.0070 



a The Ui and a s columns represent internal and systematic errors respectively. 
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Table 8. HPOL Synthetic R Filter Data 



Phase HJD - 2,400,000 %Q %U ISP sub %Q ISP sub %U (%) cr s a (%) 



Band Reticon: 












0.000 


49663 


.46 





.4242 


-0 


.5664 


0.121 


48901 


.67 


0. 


.3360 


-0. 


.4694 


0.206 


49562. 


.62 





.3777 


-0 


.3521 


0.292 


49615 


.49 





.3735 


-0 


.4713 


0.358 


49564 


.59 





3955 


-0. 


.4077 


0.368 


49603 


.53 





.3917 


-0 


.4168 


0.525 


48984 


.52 





.2297 


-0 


.4119 


0.657 


48908 


.60 


0. 


4936 


-0. 


.4751 


0.661 


48895 


.71 





.4402 


-0. 


.4586 


0.668 


49555 


.66 





.4711 


-0 


.4491 


0.737 


48922 


.58 





.3817 


-0. 


.4389 


0.762 


49194 


.60 





.4064 


-0 


.4960 


0.890 


49506 


.78 





.2308 


-0 


.5012 


0.962 


49533. 


.59 


0. 


.4678 


-0. 


.5736 


R Band CCD: 












0.000 


50685 


.64 





.4015 


-0 


.6876 


0.004 


50698 


.63 


0. 


.5123 


-0. 


.4037 


0.008 


50685 


.74 





.3668 


-0. 


.7072 


0.010 


50698 


.71 





.4717 


-0 


.3918 


0.011 


50724. 


.60 


0. 


.5048 


-0. 


.5398 


0.015 


50685 


.83 


0. 


.3437 


-0. 


.7394 


0.019 


50724. 


.70 


0. 


.4562 


-0. 


.5390 


0.039 


49935 


.68 


0. 


.4169 


-0 


.5498 


0.043 


49922 


.79 


0. 


.5921 


-0. 


.5174 


0.069 


50647 


.71 


0. 


.4405 


-0. 


.4660 


0.073 


50634 


.82 





.3695 


-0 


.3606 


0.079 


50686. 


.66 


0. 


.3101 


-0. 


.5081 


0.086 


50725. 


.57 


0. 


.3175 


-0. 


.4057 


0.087 


50712 


.64 





.3317 


-0 


.4421 


0.117 


49910 


.82 





.3510 


-0 


.3731 


0.145 


50661 


.64 





3317 


-0. 


.4227 


0.184 


49872 


.86 





.2922 


-0 


.3812 


0.219 


50636. 


.72 


0. 


.3710 


-0. 


.4600 


0.226 


50636. 


.80 


0. 


.3663 


-0. 


.4598 


0.283 


50585 


.78 





.4810 


-0 


.5997 


0.287 


51025 


.78 





.3534 


-0 


.4673 


0.300 


51064 


.76 





.3697 


-0. 


.4742 


0.333 


50922 


.86 





.4735 


-0 


.4876 


0.354 


50987 


.82 





.4001 


-0 


.4492 


0.385 


50690 


.62 


0. 


.4323 


-0. 


.4968 


0.401 


50703. 


.77 


0. 


.3302 


-0. 


.4507 


0.403 


50716 


.73 





.3535 


-0 


.4378 


0.410 


49901 


.66 





.3238 


-0. 


.2942 


0.447 


50639. 


.67 


0. 


.2888 


-0 


.3107 


0.455 


50639 


.77 





.3009 


-0 


.3000 



0.2205 


-0.2422 


0.0027 


0.0200 


0.1341 


-0.1494 


0.0038 


0.0200 


0.1861 


-0.0332 


0.0067 


0.0200 


0.1709 


-0.1498 


0.0025 


0.0200 


0.1947 


-0.0863 


0.0032 


0.0200 


0.1852 


-0.0923 


0.0027 


0.0200 


0.0295 


-0.0905 


0.0032 


0.0200 


0.2912 


-0.1619 


0.0045 


0.0200 


0.2373 


-0.1349 


0.0033 


0.0200 


0.2711 


-0.1300 


0.0022 


0.0200 


0.1567 


-0.1187 


0.0042 


0.0200 


0.2035 


-0.1761 


0.0018 


0.0200 


0.0096 


-0.1950 


0.0098 


0.0200 


0.2714 


-0.2479 


0.0035 


0.0200 


0.1993 


-0.3676 


0.0020 


0.0080 


0.3090 


-0.0828 


0.0014 


0.0080 


0.1640 


-0.3878 


0.0019 


0.0080 


0.2708 


-0.0742 


0.0019 


0.0080 


0.3034 


-0.2217 


0.0014 


0.0080 


0.1418 


-0.4221 


0.0036 


0.0080 


0.2541 


-0.2205 


0.0017 


0.0080 


0.2151 


-0.2291 


0.0014 


0.0040 


0.3832 


-0.1864 


0.0025 


0.0040 


0.2384 


-0.1469 


0.0015 


0.0070 


0.1678 


-0.0419 


0.0010 


0.0070 


0.1085 


-0.1898 


0.0011 


0.0080 


0.1147 


-0.0859 


0.0013 


0.0080 


0.1289 


-0.1217 


0.0012 


0.0080 


0.1492 


-0.0526 


0.0010 


0.0040 


0.1267 


-0.0989 


0.0015 


0.0070 


0.0882 


-0.0616 


0.0011 


0.0040 


0.1691 


-0.1402 


0.0012 


0.0070 


0.1638 


-0.1391 


0.0011 


0.0070 


0.2692 


-0.2653 


0.0015 


0.0070 


0.1512 


-0.1477 


0.0011 


0.0140 


0.1683 


-0.1560 


0.0018 


0.0140 


0.2710 


-0.1675 


0.0010 


0.0080 


0.1984 


-0.1304 


0.0010 


0.0140 


0.2300 


-0.1768 


0.0012 


0.0080 


0.1257 


-0.1271 


0.0015 


0.0080 


0.1509 


-0.1187 


0.0020 


0.0080 


0.1197 


0.0269 


0.0012 


0.0040 


0.0874 


0.0085 


0.0011 


0.0070 


0.0987 


0.0201 


0.0011 


0.0070 
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Table 8 — Continued 



Phase HJD - 2,400,000 %Q %U ISP sub %Q ISP sub %U (%) cr s a (%) 



0.463 


50639.87 


0.3125 


-0.3099 


0.1106 


0.0090 


0.0017 


0.0070 


0.483 


50316.65 


0.4170 


-0.3947 


0.2131 


-0.0749 


0.0009 


0.0200 


0.485 


49863.82 


0.2997 


-0.2970 


0.0976 


0.0219 


0.0010 


0.0040 


0.486 


50769.56 


0.3217 


-0.3832 


0.1196 


-0.0638 


0.0014 


0.0080 


0.525 


50640.68 


0.3155 


-0.4577 


0.1120 


-0.1365 


0.0014 


0.0070 


0.534 


50640.79 


0.2685 


-0.4502 


0.0701 


-0.1377 


0.0016 


0.0070 


0.559 


49864.78 


0.5560 


-0.4658 


0.3583 


-0.1510 


0.0022 


0.0040 


0.572 


49916.70 


0.4179 


-0.3970 


0.2160 


-0.0770 


0.0009 


0.0040 


0.645 


50797.50 


0.4237 


-0.4673 


0.2253 


-0.1500 


0.0015 


0.0080 


0.653 


49943.63 


0.4946 


-0.4361 


0.2744 


-0.1159 


0.0011 


0.0040 


0.671 


51056.62 


0.4657 


-0.4839 


0.2640 


-0.1650 


0.0009 


0.0140 


0.723 


50927.90 


0.4251 


-0.5203 


0.2193 


-0.1945 


0.0028 


0.0080 


0.749 


49970.74 


0.4687 


-0.4084 


0.2651 


-0.0870 


0.0009 


0.0040 


0.796 


49867.84 


0.4310 


-0.4647 


0.2274 


-0.1439 


0.0008 


0.0040 


0.833 


50644.66 


0.3140 


-0.4453 


0.1120 


-0.1262 


0.0010 


0.0070 


0.842 


50644.77 


0.2980 


-0.4499 


0.0959 


-0.1305 


0.0011 


0.0070 


0.848 


49790.89 


0.3819 


-0.4402 


0.1778 


-0.1185 


0.0014 


0.0040 


0.848 


50644.86 


0.2874 


-0.4334 


0.0851 


-0.1131 


0.0012 


0.0070 


0.862 


49842.82 


0.3860 


-0.4130 


0.1842 


-0.0936 


0.0010 


0.0040 


0.912 


50645.68 


0.3574 


-0.3987 


0.1543 


-0.0783 


0.0011 


0.0070 


0.920 


50645.78 


0.3495 


-0.3898 


0.1483 


-0.0707 


0.0014 


0.0070 


0.939 


50749.54 


0.4245 


-0.5886 


0.2219 


-0.2686 


0.0015 


0.0080 


0.955 


49908.72 


0.6552 


-0.4704 


0.4477 


-0.1407 


0.0039 


0.0040 


0.969 


49947.72 


0.5424 


-0.4346 


0.3401 


-0.1159 


0.0016 


0.0040 


0.980 


50594.80 


0.4749 


-0.5607 


0.2626 


-0.2252 


0.0059 


0.0070 



a The Ui and a s columns represent internal and systematic errors respectively. 
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Table 9. HPOL Synthetic / Filter Data 



Phase HJD - 2,400,000 %Q %U ISP sub %Q ISP sub %U of (%) cr s a (%) 



/ Band Reticon: 










0.000 


49663. 


46 


0.3626 


-0. 


5178 


0.121 


48901. 


.67 


0.3121 


-0 


.5255 


0.206 


49562. 


62 


0.2211 


-0 


3197 


0.292 


49615. 


.49 


0.2800 


-0 


.4636 


0.358 


49564. 


.59 


0.2809 


-0 


.3208 


0.368 


49603. 


53 


0.3388 


-0 


.4346 


0.525 


48984. 


52 


0.2337 


-0 


.3604 


0.657 


48908. 


.60 


0.4512 


-0 


.4725 


0.661 


48895. 


.71 


0.3869 


-0 


.4243 


0.668 


49555. 


66 


0.4397 


-0 


.3904 


0.737 


48922. 


.58 


0.2630 


-0 


.2777 


0.762 


49194. 


60 


0.3384 


-0 


.4473 


0.890 


49506. 


78 


-0.3013 


-0. 


.9340 


0.962 


49533. 


.59 


0.4211 


-0. 


.5398 


/ Band CCD: 










0.000 


50685. 


.64 


0.3804 


-0 


5442 


0.004 


50698. 


.63 


0.4507 


-0 


3146 


0.008 


50685. 


.74 


0.3583 


-0 


.5498 


0.010 


50698. 


71 


0.4507 


-0 


.3272 


0.011 


50724. 


.60 


0.4749 


-0 


.4690 


0.015 


50685. 


.83 


0.3397 


-0 


.5846 


0.019 


50724 


70 


0.4457 


-0 


4810 


0.039 


49935. 


68 


0.3865 


-0 


4621 


0.043 


49922. 


.79 


0.5194 


-0 


4285 


0.069 


50647. 


71 


0.3885 


-0 


.3872 


0.073 


50634. 


.82 


0.3375 


-0 


.3073 


0.079 


50686. 


.66 


0.2636 


-0 


.4484 


0.086 


50725. 


.57 


0.2855 


-0. 


.3658 


0.087 


50712. 


64 


0.2864 


-0 


.3721 


0.117 


49910. 


.82 


0.2913 


-0 


.3028 


0.145 


50661. 


.64 


0.2625 


-0 


3677 


0.184 


49872. 


86 


0.2252 


-0 


.2967 


0.219 


50636. 


.72 


0.2874 


-0 


.3558 


0.226 


50636. 


.80 


0.2898 


-0 


.3749 


0.287 


51025. 


.78 


0.2746 


-0 


3784 


0.300 


51064 


76 


0.2952 


-0 


.4062 


0.333 


50922. 


.86 


0.3977 


-0 


.4008 


0.354 


50987. 


.82 


0.3203 


-0. 


.3821 


0.385 


50690. 


62 


0.3443 


-0 


.3991 


0.401 


50703. 


.77 


0.2546 


-0. 


.3997 


0.403 


50716. 


.73 


0.2897 


-0 


.3675 


0.410 


49901. 


66 


0.2705 


-0 


.2435 


0.447 


50639. 


.67 


0.2618 


-0. 


.2654 


0.455 


50639. 


.77 


0.2673 


-0 


.2538 


0.463 


50639. 


87 


0.2749 


-0 


.2538 



0.1859 


-0.1806 


0.0107 


0.0200 


0.1643 


-0.2718 


0.0210 


0.0200 


0.2852 


-0.0712 


0.0212 


0.0200 


0.1333 


-0.1712 


0.0147 


0.0200 


0.1720 


-0.0486 


0.0135 


0.0200 


0.1126 


-0.0753 


0.0167 


0.0200 


0.1426 


-0.1002 


0.0127 


0.0200 


0.2977 


-0.3356 


0.0278 


0.0200 


0.2265 


-0.0991 


0.0178 


0.0200 


0.3286 


-0.1580 


0.0090 


0.0200 


-0.2879 


-0.0631 


0.0139 


0.0200 


0.1821 


-0.1987 


0.0061 


0.0200 


-0.4458 


-0.6993 


0.0649 


0.0200 


0.3808 


-0.1714 


0.0138 


0.0200 


0.2091 


-0.2734 


0.0019 


0.0080 


0.2795 


-0.0439 


0.0014 


0.0080 


0.1870 


-0.2790 


0.0018 


0.0080 


0.2794 


-0.0564 


0.0018 


0.0080 


0.3036 


-0.1983 


0.0014 


0.0080 


0.1683 


-0.3138 


0.0038 


0.0080 


0.2744 


-0.2102 


0.0017 


0.0080 


0.2152 


-0.1913 


0.0015 


0.0070 


0.3481 


-0.1578 


0.0016 


0.0070 


0.2171 


-0.1163 


0.0012 


0.0080 


0.1662 


-0.0366 


0.0011 


0.0080 


0.0923 


-0.1777 


0.0012 


0.0080 


0.1142 


-0.0951 


0.0013 


0.0080 


0.1150 


-0.1013 


0.0012 


0.0080 


0.1198 


-0.0321 


0.0010 


0.0070 


0.0912 


-0.0968 


0.0019 


0.0080 


0.0538 


-0.0259 


0.0012 


0.0070 


0.1159 


-0.0849 


0.0012 


0.0080 


0.1183 


-0.1040 


0.0011 


0.0080 


0.1034 


-0.1078 


0.0011 


0.0220 


0.1240 


-0.1356 


0.0017 


0.0220 


0.2263 


-0.1298 


0.0011 


0.0060 


0.1489 


-0.1112 


0.0010 


0.0220 


0.1731 


-0.1284 


0.0011 


0.0080 


0.0833 


-0.1289 


0.0013 


0.0080 


0.1184 


-0.0967 


0.0020 


0.0080 


0.0991 


0.0274 


0.0013 


0.0070 


0.0902 


0.0056 


0.0012 


0.0080 


0.0959 


0.0171 


0.0012 


0.0080 


0.1035 


0.0171 


0.0013 


0.0080 
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Table 9 — Continued 



Phase 


HJD - 2,400,000 


%Q 


%U 


ISP sub %Q 


ISP sub % U 


<T S a (%) 


CTs a (%) 


0.483 


50316.65 


0.3544 


-0.3245 


0.1829 


-0.0535 


0.0010 


0.0200 


0.485 


49863.82 


0.2575 


-0.2304 


0.0861 


0.0405 


0.0012 


0.0070 


0.486 


50769.56 


0.2740 


-0.3214 


0.1030 


-0.0509 


0.0014 


0.0080 


0.525 


50640.68 


0.3033 


-0.3901 


0.1319 


-0.1192 


0.0018 


0.0080 


0.534 


50640.79 


0.2402 


-0.3920 


0.0689 


-0.1213 


0.0013 


0.0080 


0.559 


49864.78 


0.5150 


-0.4246 


0.3438 


-0.1540 


0.0022 


0.0070 


0.572 


49916.70 


0.3529 


-0.3019 


0.1815 


-0.0309 


0.0011 


0.0070 


0.645 


50797.50 


0.3609 


-0.3924 


0.1897 


-0.1218 


0.0016 


0.0080 


0.653 


49943.63 


0.3998 


-0.3541 


0.2285 


-0.0833 


0.0011 


0.0070 


0.671 


51056.62 


0.3735 


-0.4160 


0.2021 


-0.1452 


0.0010 


0.0220 


0.723 


50927.90 


0.2960 


-0.4475 


0.1246 


-0.1766 


0.0025 


0.0060 


0.749 


49970.74 


0.3880 


-0.3424 


0.2166 


-0.0715 


0.0009 


0.0070 


0.796 


49867.84 


0.3552 


-0.3721 


0.1839 


-0.1013 


0.0009 


0.0070 


0.833 


50644.66 


0.2470 


-0.3595 


0.0756 


-0.0886 


0.0011 


0.0080 


0.842 


50644.77 


0.2378 


-0.3646 


0.0665 


-0.0937 


0.0012 


0.0080 


0.848 


49790.89 


0.2912 


-0.3398 


0.1198 


-0.0689 


0.0018 


0.0070 


0.848 


50644.86 


0.2267 


-0.3582 


0.0553 


-0.0873 


0.0014 


0.0080 


0.862 


49842.82 


0.3112 


-0.3377 


0.1400 


-0.0670 


0.0011 


0.0070 


0.912 


50645.68 


0.3068 


-0.3286 


0.1354 


-0.0577 


0.0014 


0.0080 


0.920 


50645.78 


0.2922 


-0.3419 


0.1210 


-0.0712 


0.0012 


0.0080 


0.939 


50749.54 


0.3517 


-0.4706 


0.1803 


-0.1997 


0.0014 


0.0080 


0.955 


49908.72 


0.5040 


-0.3810 


0.3327 


-0.1102 


0.0059 


0.0070 


0.969 


49947.72 


0.4425 


-0.3480 


0.2715 


-0.0772 


0.0015 


0.0070 



a The <Ji and a s columns represent internal and systematic errors respectively. 
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Table 10. Position Angle Information for /3 Lyrae 



Band 


HPOL Reticon (deg) 


HPOL CCD (deg) 


AH (deg) 


FCO (deg) 


B 


139.8 ± 3.5 


165.0 ± 0.7 


169.8 ± 1.7 


165.1 ± 0.9 


V 


147.3 ± 3.0 


164.0 ± 0.5 


171.8 ± 1.0 


163.1 ± 0.4 


R 


158.4 ± 0.4 


158.3 ± 0.4 




162.2 ± 0.7 


I 


57.3 ± 3.0 


162.7 ± 0.6 






Index 


158.8 ± 0.6 


142.1 ± 1.7 






Ha 




47.9 ± 9.7 






H/3 




71.1 ± 2.5 






He IA5876 




63.8 ± 6.4 






He IA7065 




51.0 ± 3.2 






He IA6678 




138.2 ± 6.6 







Note. - - These values were calculated by fitting a line in Q-U space to the data, 
excluding points between phases 0.475 and 0.525. For the HPOL data, we used system- 
atic errors in both Q and U in the least-squares fit calculation for all bands except the 
Balmer jump index (the vector difference between the polarization above and below 
the Balmer jump) and the lines, for which we used the intrinsic errors. See Sections 
3.1 and 3.2. 
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Table 11. BVRI Fourier Fit Parameters 



Band 


Zero Point 


Term 


n 


A 




B 


0.239470 


1 


2.98 ± 0.01 


0.127 ± 0.009 


0.30 ± 0.01 






2 


1.00 ± 0.02 


0.078 ± 0.008 


0.37 ± 0.02 


V 


0.243419 


1 


2.98 ± 0.01 


0.096 ± 0.007 


0.30 ± 0.01 






2 


1.01 ± 0.03 


0.050 ± 0.006 


0.36 ± 0.02 


R 


0.241148 


1 


3.00 ± 0.01 


0.107 ± 0.008 


0.32 ± 0.01 






2 


4.02 ± 0.03 


0.052 ± 0.009 


0.08 ± 0.03 


I 


0.195021 


1 


3.00 ± 0.01 


0.093 ± 0.006 


0.30 ± 0.01 






2 


4.01 ± 0.02 


0.050 ± 0.006 


0.10 ± 0.02 



Note. - These parameters describe our Fourier fits to the broadband 
polarimetric data in the combined HPOL, FCO and AH data sets. See 
Section 3.1. 



Table 12. Broadband Secondary Eclipse Phases and Maximum Hot Spot Size 



Band 


Phase 


Min PA Phase Rot 


Max PA Phase Rot 


Hot Spot Size in Phase 


if DC 






HSpa 


HSsm (i = 90°) 


HSsm (i = 86°) 


B 


0.481 ± 0.004 


0.448 


0.543 


0.095 


0.186 


0.045 


22 Rq 


33 R e 


31 Rq 


26 Rq 


V 


0.481 ± 0.004 


0.385 


0.554 


0.169 


0.184 


0.080 


22 Rq 


58 Rq 


18 Rq 


9 Rq 


R 


0.449 ± 0.005 


0.448 


0.524 


0.076 


0.171 


0.080 


33 Rq 


26 Rq 


19 Rq 


10 Rq 


I 


0.449 ± 0.004 


0.448 


0.524 


0.076 


0.140 


0.077 


33 Rq 


26 R e 


12 Rq 


2 Rq 



Note. — Columns two through six are parameters used in the three hot spot size estimates (see Sections 4.1, 4.2, and 4.3). Flux normalized to maximum light was used for 
if DC an d if min- The last four columns represent the width of the hot spot calculated using our different estimates (see Section 4). All three polarimetric data sets were used 
in this analysis and the resulting estimates. 
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Fig. 1. — Data points represent the B band polarimetric observations from HPOL Reticon 
(diamonds), HPOL CCD (squar es), AH (circles) and FC O (triangles). From top: Normalized 
V band Fourier fit light curve (IHarmanec et al.lll996l ). projected polarization (see Section 
3.1), and position angle (degrees) versus phase. Error bars are shown for uncertainties larger 
than 0.025 in %Q p and 5.0° in position angle. The HPOL error bars shown represent the 
larger of the intrinsic and systematic uncertainties. All data have been wrapped so that more 
than one complete period is shown. The solid line in the middle panel represents our Fourier 
fit to the %Q p data (see Section 3.1). The dotted line represents zero projected polarization. 
The solid line in the bottom panel represents the average position angle of 164°. 
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Fig. 2. — Same as Figure 1, but for V band polarimetry. The dashed blue curve (in the 
online version, otherwise dashed light grey) represents the three-term Fourier fit and the red 
curve (in the online version, otherwise solid) represents the two-term Fourier fit. Note the 
two fits differ at primary and secondary eclipse, the quadrature phases, and near phases 0.15 
and 0.9. 
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Fig. 3. — Same as Figure 1, but for R band polarimetry. 
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Fig. 4. — Same as Figure 1, but for / band polarimetry. 
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Fig. 5. — From top: Normalized V band Fourier fit light curve and projected polarized flux 
curves for the BVRI bands (see Section 3.1). The projected polarized flux curves are formed 
by multiplying each band's Fourier fit polarization curve by its normalized to maximum light 
Fourier fit light curve. 
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He I A5876 emission line polarization from HPOL C CD observations (see S ection 

1996h . pro- 



Fig. 6. 

3.2). From top: Normalized V band Fourier fit light curve (IHarmanec et al. 



jected polarization, and position angle (degrees) versus phase. Error bars represent intrinsic 
uncertainties. All data have been wrapped in phase to display more than one complete 
period. 
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Fig. 7. — Proposed geometry of the (3 Lyr system at various phases in our proposed hot spot 
model. Arrows represent the polarization arising from the disk edge. From Top: The first 
maximum in polarization occurs at the first quadrature phase. The polarization then begins 
to decrease as the hot spot (hatched region) rotates into view. The minimum in polarization 
occurs when the area eclipsed by the loser and disrupted by the hot spot is maximized. The 
minimum in flux occurs as the hot spot is rotating off the edge of the visible disk. The 
second maximum in polarization occurs at the second quadrature phase. 
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Fig. 8. — Sketch of the geometry we used for our hot spot size estimates (not to scale, 
see Section 4). We used parameters obtained by Linnell (2000) for the disk height, disk 
diameter, loser radius and separation between the two components. The hatched region in 
the Observer's View represents the hot spot size based on the %Q p method (see Section 4.1). 
The blackened area of the circle in the Observer's View is the uneclipsed area of the loser at 
primary eclipse. The filled square inside the accretion disk in the Polar View indicates the 
location of the roots of the bipolar outflows as given by the H a absorption core whil e the 
filled star represents the same thing for the Ha emission wings (IHarmanec et al.lll996l ). 
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Fig. 9. — Sketch of the geometry we used for our simple model hot spot size estimate at an 
inclination angle of 86° (not to scale; see Section 4.3). The hatched regions represent the 
additional area of the disk in our estimate due to the assumption of a rectangularly shaped 
disk. 



